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Abstract 


Flat-plate  penetrometers,  such  as  the  Marchetti  Dilatometer  and  K0  Stepped- Blade, 
have  become  an  important  part  of  in  situ  testing  in  geotechnical  engineering.  However, 
the  use  of  flat-plate  penetrometers  has  been  highly  empirical,  mainly  due  to  the  lack 
of  knowledge  of  soil  response  to  the  flat-plate  penetration.  In  this  research  project, 
a  numerical  technique  capable  of  computing  strain  paths  for  three-dimensional  pen¬ 
etrometers  was  developed.  A  calibration  chamber  system  for  cohesive  soils  and  model 
flat-plates  were  fabricated.  For  the  first  time,  three-dimensional  strain  path  analyses 
were  performed  for  several  of  the  flat-plate  penetrometers  currently  being  used  in 
the  industry.  The  results  show  that  flat-plates  can  induce  large  strains  and  strain 
reversals  at  levels  comparable  to  those  of  cone  penetration.  The  characteristics  of  the 
strain  field  during  a  flat-plate  penetration  is  influenced  by  both  the  w/t  ratio  and  the 
tip  apex  angle.  Regardless  of  the  geometry,  the  pore  pressure  and  total  stress  peak  at 
the  tip  of  the  penetrometer.  The  pore  pressure  and  total  stress  decrease  sharply  as  the 
soil  element  passes  the  tip  of  the  flat-plate.  The  pore  presure  and  total  stress  around 
the  flat-plate  do  not  increase  with  the  plate  thickness  as  some  had  suggested.  For 
a  simple  flat-plate  (i.e.,  the  flat  Dilatometer),  the  penetration  induced  pore  pressure 
is  positively  related  to  the  soil  rigidity  index  in  the  horizontal  direction.  The  excess 
pore  pressure  dissipation  around  the  flat-plate  in  an  anisotropically  consolidated  clay 
follows  an  axisymmetric  pattern.  Interpretation  of  the  pore  pressure  dissipation  data 
using  an  axisymmetric  solution  resulted  a  C\  value  comparable  but  lower  than  those 
obtained  from  reference  consolidation  tests. 
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Chapter  1 


INTRODUCTION 


1.1  Applications  of  Flat-Plate  Penetrometers  in 
Geotechnical  Engineering 

The  technique  of  inserting  a  flat-pla4e  into  soil  to  obtain  engineering  design  param¬ 
eters  has  been  incorporated  in  a  number  of  recent  insitu  soil  testing  techniques. 
These  include  the  Marchetti  (1975;  1980)  Dilatometer  (DMT),  the  AVStepped  Blade 
(Handy  et  al.,  1982;  Lutenegger  and  Timian,  1986),  and  various  push-in  spade  cells 
(Massarsch,  1975;  Tedd  and  Charles,  1981).  Figure  1.1  presents  a  general  description 
of  some  of  the  commonly  used  flat-plates.  Flat-plate  penetrometers  offer  potentially 
important  advantages  over  other  (i.e.  cone  penetrometer)  penetrometers  including: 


1.  It  is  possible  to  reveal  directional  anisotropy  in  horizontal  stress  (Handy  et  al., 
1982); 

2.  The  penetration  simulates  boundary  conditions  in  the  case  of  a  one-dimensional 
compression,  therefore  enables  a  simplified  interpretation  of  the  test  data  (Marchetti, 
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Steppe!  Ht4c  Ttferti  lUle 


lilitaeter 


Spile  Cells 


Figure  1.1:  Some  of  the  commonly  used  flat-plate  penetrometers  (Handy  et  al.  1982; 
Marchetti,  1980;  Handy  and  Lutenegger,  1985;  Ladd  et  al.  1979;  and  Fukuoka  and 
Imamura,  1983)  2 


1980;  Handy  et  al.,  1982); 


3.  A  sharp  and  relatively  thin  flat-plate  penetrometer  causes  less  soil  disturbance 
(as  opposed  to  a  60  degree  cone  penetrometer)  and  therefore  requires  less  ex¬ 
trapolation  in  predicting  the  undisturbed  soil  parameters  (Marchetti,  1980); 

4.  The  tests  are  relatively  simple  to  perform,  cost  effective  and  results  are  poten¬ 
tially  useful  (Schmertmann,  1981). 

One  might  also  add  that  it  is  mechanically  easier  to  install  pressure  sensors  on  the 
face  of  a  flat-plate  as  evidenced  by  the  designs  of  the  Dilatometer  and  the  /^-Stepped 
Blade. 

In  a  very  general  sense,  the  flat-plate  penetration  can  be  considered  a  special  class 
of  soil  failure  problems,  i.  e.  ,  the  soil  is  failed  during  insertion  to  give  room  for  the 
penetrometer.  Scott  (1987)  has  suggested  that  there  are  three  components  of  a  failure 
and  for  it  to  be  understood  in  the  most  general  sense,  information  is  necessary  on 
all  three  components.  They  are  mechanism,  properties,  and  analysis.  A  common 
approach  in  interpreting  flat-plate  penetration  tests  is  to  assume  the  mechanisms  (e. 
g.  ,  plane  strain  or  cylindrical  cavity  expansion),  perform  the  analysis  based  on  the 
test  data  (e.  g.  ,  tip  resistance  or  membrane  expansion  pressure)  and  back  calculate 
the  soil  properties  (e.  g.  ,  shear  strength  or  modulus).  The  validity  of  this  approach 
is  directly  related  to  that  of  the  assumed  mechanisms  which  are,  in  general,  the  least 
known. 


1.2  Lateral  Compression  in  An  Elastic  Half  Space 

Handy  et  al.  ,  (1982)  admitted  the  lack  of  knowledge  regarding  the  soil  disturbance 
during  penetration  and  stated  that  ”. . .  Although  disturbance  cannot  be  directly  eval¬ 
uated,  it  should  relate  to  thickness  of  a  penetrating  blade  . . . ".  The  AVStepped  Blade 
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was  actually  conceived  based  on  this  premise.  Similar  conjectures  were  used  to  form 
the  basis  for  many  other  types  of  flat-plate  penetrometers.  The  situation  is  under¬ 
standable  because  of  the  large  strain  and  three-dimensional  nature  of  flat-plate  pen¬ 
etrations.  Essentially,  all  the  currently  available  interpretation  methods  (Marchetti, 
1980;  Handy  et  al.  ,1982;  Tedd  and  Charles,  1983)  for  flat-  plate  penetration  tests  are 
based  on  elastic  theories.  They  deal  with  a  boundary  value  problem  of  a  rigid  plate 
translating  laterally  in  an  elastic  half  space  (Figure  1.2).  Liam  Finn  (1963)  solved 
the  problem  considering  the  plate  as  infinitely  long.  The  stress  at  the  center  of  the 
rigid  plate  is  related  to  displacement  as: 


d  E 

’v  "  5^(1  -i/2) 


(1.1) 


This  is  the  basis  for  the  interpretation  of  spade  cell  tests  a s  suggested  by  Tedd  and 
Charles  (1983).  Marchetti  (1980)  used  a  solution  essentially  for  the  same  boundary 
conditions  except  the  plate  is  circular  and 

rd  E 


Py  2D  (1  —  v2) 

For  the  DMT,  D  =  60  mm,  d  =  1  mm,  Eq.  1.2  becomes: 

E 


=  38.2  (P‘ 


DMT 


-P, 


DMT 


) 


(1.2) 


(1.3) 


The  left  hand  term  of  Eq.  1.3  is  commonly  referred  to  as  the  dilatometer  modulus, 
Ed. 


Handy  et  al.  (1982)  adopted  a  suggestion  by  Janbu  (1967)  that  the  soil  compressive 
modulus  is  proportional  to  stress.  An  exponential  relationship  between  Py  and  d  (= 
t/2)  was  derived  as: 

Py  =  P0aebt  (1.4) 


Eq.  1.4  provides  the  basis  of  extrapolating  the  insitu  lateral  earth  pressure  from  A'o- 
Stepped  Blade  tests. 
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Figure  1.2:  Lateral  compression  in  elastic  half  space  (Liam  Finn,  1963) 
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Although  never  clearly  stated  before,  the  use  of  lateral  compression  approach  in  an¬ 
alyzing  the  flat-plate  penetration  problems  is  obviously  influenced  by  the  cylindrical 
cavity  expansion  theory  (Vesic,  1972).  It  considers  the  penetration  mechanisms  as 
such  that  strains  in  the  horizontal  direction  are  the  predominant  component  during 
a  vertical  penetration.  For  cylindrical  penetrometers,  i.e.  ,  piles  or  cones,  Baligh 
(1985)  had  shown  that  the  lateral  strains  do  represent  the  predominant  component. 
However,  there  has  been  little  evidence  to  prove  that  the  same  is  true  for  flat-plate 
penetrations.  Had  this  not  been  the  case,  then  the  lateral  compression  approach  is  not 
even  relevant  to  the  problem  it  intends  to  solve  except  for  the  case  of  DMT  expansion. 

A  fundamental  understanding  of  how  soils  respond  to  the  insertion  of  a  flat-plate 
penetrometer  is  paramount  to  the  success  of  interpreting  the  results  produced  by 
such  instruments  and  the  continued  application  to  Geotechnical  problems.  Without 
a  rational  basis,  the  interpretation  of  results  will  remain  strictly  empirical  and  will 
be  limited  by  local  experience  with  local  soils. 


1.3  Objectives  of  The  Research 

The  objectives  of  this  research  were  to  : 

1.  Estimate  the  soil  strain  field  during  a  steady  penetration  of  a  flat-plate  insertion 
in  saturated  cohesive  soils  and  study  its  relationship  with  the  geometry  of  a 
penetrometer; 

2.  Evaluate  the  strain  components  and  unique  features  of  the  strain  field  in  the 
soil  during  the  flat-plate  penetration; 
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3.  Perform  a  general  study  of  the  total  stress  and  pore  pressure  distribution  in  the 
soil  as  a  result  of  a  flat-plate  penetration; 

4.  Validate  the  results  of  the  above  studies  by  performing  laboratory  flat-plate 
penetration  tests  in  a  calibration  chamber; 

5.  Conduct  laboratory  experiments  on  the  same  soil  used  in  the  chamber  testing 
to  provide  reference  soil  properties.  Detailed  studies  regarding  the  soil  response 
to  flat-plate  induced  strain  paths  will  then  follow; 

6.  Evaluate  available  field  test  data  in  light  of  the  analytical  and  experimental 
studies. 

To  achieve  the  objectives,  the  research  was  divided  into  three  tasks.  The  first  task  con¬ 
centrated  on  analytical  studies.  A  numerical  technique  which  is  capable  of  performing 
strain  path  analyses  for  arbitrary  three-  dimensional  penetrometers  was  developed. 
The  corresponding  stress  field  was  determined  using  a  Modified  Cam  Clay  Model. 
The  pore  pressure  distribution  was  estimated  by  considering  the  stress  equilibrium. 
A  series  of  reference  laboratory  experiments  and  penetration  tests  of  a  model  flat 
Dilatometer  were  performed  in  the  second  task.  A  calibration  chamber  system  along 
with  other  testing  equipment  associated  with  the  research  were  fabricated.  The  cali¬ 
bration  chamber  was  used  to  perform  model  flat  Dilatometer  penetration  tests  under 
a  controlled  environment. 

A  series  of  field  flat-plate  penetration  test  data  along  with  other  in  situ  and  labora¬ 
tory  test  results  from  a  research  test  site  in  Massena,  New  York  were  compiled  and 
studied  in  the  third  task.  The  report  describes  the  details  of  the  research,  presents 
the  results  and  discusses  the  merits  of  flat-plate  penetration  tests  as  a  result  of  the 
research. 
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Chapter  2 


Three-Dimensional  Strain  Path 
Analyses 


2.1  Background 

The  strain  path  method  (SPM)  was  developed  by  M.M.  Baligh  (1985)  as  a  result  of 
his  research  in  deep  foundations.  The  method  was  inspired  by  the  observations  of 
soil  deformations  caused  by  the  undrained  penetration  of  rigid  objects  in  saturated 
clays.  Baligh  (1985)  hypothesized  that,  ’’...due  to  the  severe  kinematic  constraints 
that  exist  in  ’’deep”  penetration  problems,  soil  deformations  and  strains  are,  by  and 
large,  independent  of  the  shearing  resistance  of  the  soil.  This  means  that  these  prob¬ 
lems  are  essentially  strain-controlled  and  implies  that,  even  if  relatively  "simple”  soil 
properties  (e.  g.  ,  isotropy)  are  utilized  to  estimate  deformations  and  strains  caused 
by  penetration,  the  errors  introduced  are  expected  to  be  small.  ...”.  Based  on 
these  premises  and  further  assumptions  that  a  saturated  clay  is  inviscid  and  incom¬ 
pressible,  the  strains  are  decoupled  from  stresses;  a  result  that  simplifies  the  strain 
path  analysis  to  a  potential  flow  problem.  Valuable  contributions  have  been  made 
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by  Baligh  and  his  coworkers  using  SPM  in  the  understanding  of  the  soil  response 
due  to  penetration  of  solid  as  well  as  hollow,  cylindrical  penetrometers  in  saturated 
clays  (Baligh,  1985;  1986;  1986a;  and  Baligh,  et  al.  ,  1987).  A  similar  technique 
was  also  developed  by  Tumay,  et  al.  (1985)  using  conformal  mapping.  The  SPM  is 
superior  to  the  conventional  cavity  expansion  theories  (Vesic,  1972)  in  many  ways 
as  indicated  by  Baligh’s  work.  The  current  SPM  is,  however,  applicable  only  to  ax- 
isymmetric  penetrometers  such  as  a  cone  or  a  tube  sampler.  Calculations  of  potential 
flow  about  an  axisymmetric  or  essentially  any  two-dimensional  bodies  are  well  estab¬ 
lished  in  fluid  dynamics  and  solutions  are  readily  available  (Chow,  1982).  To  apply 
the  SPM  concept  to  flat-plate  penetrations,  however,  the  above  mentioned  techniques 
have  to  be  extended  to  three-dimensional  cases.  A  literature  review  indicated  that 
such  a  challenge  was  not  unique  to  Geotechnical  Engineers.  The  need  to  calculate 
fluid  flows  about  bodies  of  complicated  configurations  such  as  an  airplane  with  wings, 
tail,  and  engine  appendages  stimulated  some  aeronautical  engineers  and  developed 
a  numerical  technique  commonly  known  as  the  "panel  method”  for  handling  bodies 
of  arbitrary  shape  (Chow,  1982;  Hess  and  Smith,  1962;  1964;  and  Hess,  1974).  A 
numerical  technique  was  developed  based  on  this  method  in  the  three-dimensional 
strain  path  analyses  for  flat-plate  penetrations.  The  following  is  a  description  of  the 
theoretical  background  and  development  of  this  technique. 


2.2  Mathematical  Statement  of  Potential  Flow 
Problem 

The  problem  considered  is  that  of  the  steady  flow  of  an  ideal  fluid  (inviscid  and 
incompressible)  about  a  three-dimensional  body.  Consider  the  surface  of  this  body, 
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S  be  described  as 


F(xi)  =  0  (2.1) 

where  a-,  are  Cartesian  coordinates.  The  onset  flow,  i.  e.  ,  the  flow  in  the  fluid  before 
disturbance  by  the  body,  is  taken  as  a  uniform  stream  of  unit  magnitude,  and  is 
represented  as  Vi  and 

I  VT  I  (2-2) 

The  condition  of  a  steady,  onset  flow  and  the  intrusion  of  a  three-dimensional  body 
is  shown  in  Figure  2.1.  In  the  case  of  SPM,  the  onset  flow  velocity  is  equivalent 
to  the  steady  penetration  rate  of  a  penetrometer.  The  fluid  velocity  K  at  a  point 
may  be  expressed  as  the  negative  gradient  of  a  potential  function  <I>,  which  has  two 
components 

V,  =  =  -dt(p°°  +  p)  (2.3) 

where 

P°°  =  -  /  vrdxi  (2.4) 

is  the  uniform  stream  potential  associated  with  the  onset  flow. 

P  =  -jvrdxl 

or 

Vi  =  -dip  (2.5) 

p  is  the  disturbance  potential  and  Vi  is  the  induced  flow  velocity  due  to  the  body. 
The  $  satisfies  the  Laplace’s  equation  in  the  region  R  that  is  exterior  to  S,  has  a  zero 
normal  derivative  on  S  (Figure  2.1),  and  approaches  the  uniform  stream  potential  at 
infinity.  Symbolically 

did#  =  dtdt(p°°  +  p)  =  0  (2.6) 

Since  the  onset  flow  is  uniform 

dAp°°  =  0 
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so  from  Eq.  2.6 


d,dip  =  0 


The  fact  that  the  body  surface  is  impervious  requires  that 


riidi$  |f=o=  0  on  S 


(2.7) 


(2.8) 


or 


nAp  |f=o=  -niV° 


lim  $  =  p° 

X{X%—*00 


(2.9) 


or 


lim  p  =  0 

x,r,— »oo 

since  nt  is  the  unit  outward  normal  vector  at  a  point  on  S  (Figure  2.1) 

d,F 


n,  = 


d,F 


(2.10) 


JF=0 


If  the  body  surface  is  represented  by  a  series  of  sources  q  (Figure  2.1),  each  having  a 
density  $(<7),  then  the  induced  potential  p  at  a  point  P(xi)  is  (Kellog,  1929) 


*■> 


(2.11) 


where  r(P,q)  is  the  distance  between  points  P  and  q.  The  p  as  given  in  EQ.  2.11 
satisfies  EQ.  2.6  and  EQ.  2.9.  The  condition  of  equation  EQ.  2.8  is  used  to  solve 
for 

-  nidi$  |F=0=  2*c{q)  -  /  /  +  n,-^00  =  0  (2.12) 

J  Js  r(g,q ) 


or 


2 7ro(q)  -(j  dS  =  -nxvr  (2.13) 

J  Js  r(g,q ) 

In  EQ.  2.12,  g  represents  a  point  on  the  body  surface  S.  The  second  term  on  the 
left  hand  side  of  EQ.  2.13  reaches  infinity  as  r(g,q )  — ►  0  or  as  g  approaches  q.  The 
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numerical  solution  as  described  later  ensures  that  as  g  approaches  q,  EQ.  2.13  reduces 
to 

2na(q)  +  mV°°  =  0  (2.14) 

EQ.  2.13  is  a  two-dimensional  Fredholm  integral  equation  of  the  second  kind.  Once 
<?(q)  are  solved  from  EQ.  2.13,  the  potential  function  p  can  be  derived  from  EQ. 
2.11  and  hence  the  induced  velocities  are  determined  following  EQ.  2.5. 

2.3  The  Panel  Method 

The  panel  method  is  essentially  a  numerical  technique  to  obtain  an  approximate 
solution  of  EQ.  2.13.  The  following  sections  provide  a  brief  description  of  this  method. 
Readers  are  referred  to  Hess  and  Smith  (1962  and  1964)  for  further  details. 

2.3.1  Approximation  of  the  body  surface 

The  body  surface  is  represented  by  a  large  number  of  plane  quadrilateral  elements 
("panels”)  as  shown  in  Figure  2.2  ,  within  each  of  which  the  source  density  a  is 
assumed  constant.  The  purpose  is  to  eventually  replace  the  integral  EQ.  2.13  by 
a  set  of  linear  algebraic  equations.  The  accuracy  of  the  calculation  depends  on  the 
number  of  elements  and  the  manner  in  which  these  elements  are  distributed  over  the 
body  surface.  Generally  speaking,  elements  should  be  concentrated  in  regions  of  high 
curvature  of  the  body  surface. 
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Figure  2.2:  Approximation  of  body  surface. 
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2.3.2  Induced  flow  velocity  due  to  a  plane  element 

Following  equation  EQ.  2.11  the  induced  potential  at  a  point  P(x:)  due  to  a  single 
plane  quadrilateral  element  with  a  uniform  source  density  is 

p(x')=aIL^)  (2l5) 

and  the  induced  velocity 

VM.-M't.'j  (2.16) 

where  =  coordinates  of  point  q  on  the  plane  quadrilateral. 


To  facilitate  numerical  computation  of  EQ.  2.16,  a  local  coordinate  system  is  used  so 
that  the  quadrilateral  lies  in  the  ii,x2-coordinate  plane  as  shown  in  Figure  2.3.  Let 
the  local  coordinates  at  the  four  corners  of  the  quadrilateral  be  ct(k )  (k  =  1,2, 3, 4). 
The  corners  are  numbered  in  increasing  order  as  the  perimeter  of  the  quadrilateral  is 
traversed  in  the  clockwise  direction  as  seen  from  the  unit  normal  vector  to  the  plane 
of  quadrilateral.  Considering  a  unit  source  density  (a  =  1),  at  any  point  P{xt)  in 
the  local  coordinate  system,  the  following  induced  velocity  components  I,  are  derived 
(Hess  and  Smith,  1964) 


/  _  V-  ~  I  [r(0  ±  r(-0  ~ 

it  l  J  [r(t)  +  rO)  +  d(t,;') 

,  _  A  ci(t')  -  Cj (j )  r(i)  +  r(j)  -  d(ij) 

2  it  <*(*>»  ”  1/(0  +  r(i)  +  d{t,j) 

mn(i,j)e(i)  -  hji)\  _  ( ^(O-OO-O  ~  Hj)\ 

/  a  V  x3r{j)  ). 


where 


(2.17) 

(2.18) 
(2.19) 


d(iJ)  =  [(ci(;)  -  c,(0)2  +  (c2{j)  —  c2(*))2] 1/2  (2.20) 
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Figure  2.3:  A  plane  quadrilateral  lying  in  xi,x2-plane. 
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...  ./•  ■  \  c2(j)-c2(i) 

mn(*>  J )  /  -  \  / 

Cl(j)-C!{l) 

(2.21) 

Xi  -  Ci(t'))2  +  (x2  -  c2(t))2  +  x3]1/2 

(2.22) 

e(i)  =  x2  +  (x!  -  c,(t))2 

(2.23) 

h(i)  =  (x2  -  c2(i))(xj  -  c,(t)) 

(2.24) 

(2.25) 

and  j  =  i  +  1,  if  i  =  4  then  j  =  1. 

EQ.  2.17  to  EQ.  2.19  are  valid  essentially  for  any  polygon  plane  elements  (Hess  and 
Smith,  1964)  such  as  triangles. 

2.3.3  The  null  point 

On  each  quadrilateral  element  a  point  is  selected  where  the  flow  velocity  normal  to  the 
element  is  required  to  vanish  and  where  the  source  density  is  eventually  calculated. 
The  point,  designated  as  the  null  point,  is  taken  where  the  element  gives  to  no  velocity 
in  its  own  plane.  Analytically,  it  means  that  if  the  condition  of  EQ.  2.8  is  applied  at 
a  null  point,  no  singularity  will  occur  to  the  second  term  of  EQ.  2.13.  In  the  local 
coordinate  system,  x3  =  0  and  hence  I3  =  0  (EQ.  2.19)  for  the  null  point  since  it  lies 
on  the  element  plane.  The  determination  of  the  null  point  is  essentially  to  select  a 
local  coordinate  (xi,x2,0)  so  that  I\  —  /2  =  0  or  to  solve  the  nonlinear  simultaneous 
equations  of  EQ.  2.17  and  EQ.  2.18. 

2.3.4  Formulation  of  linear  equations 

Selection  of  null  points  enables  the  evaluation  of  EQ.  2.13  at  a  finite  number  of 
locations  on  the  approximated  body  surface.  Before  establishing  the  linear  equations, 
a  global  coordinate  system  is  needed  where  all  the  induced  velocity  components  are 
eventually  summed  and  evaluated.  Let  Vi(j,k)  be  the  velocity  vector,  in  the  global 
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coordinate,  at  the  null  point  of  jth  element  due  to  a  unit  source  density  on  the  kth 
element.  This  is  done  by  first  converting  the  jth  null  point  to  the  local  coordinate  of 
kth  element.  The  induced  velocities  I{  are  calculated  according  to  EQ.  2.17  to  EQ. 
2.19.  A  coordinate  transformation  is  then  performed  to  convert  the  induced  velocity 
into  the  global  coordinate. 

V,  =  tijlj  (2.26) 

where  £ij  is  the  direction  cosines  of  the  local  coordinates. 

The  normal  velocity  component  A(j,k)  at  the  null  point  of  jth  element  by  a  unit 
source  density  on  the  kth  element  is  the  dot  product  of  Vi(j,k)  with  the  unit  normal 
vector  to  the  jth  element  in  global  coordinates,  i.  e.  , 

A(iJ)  =  n,(j)Vt(j,k)  (2.27) 

The  normal  velocity  at  the  jth  null  point  due  to  the  entire  approximated  body  surface 
must  equal  the  negative  of  the  normal  component  of  the  onset  flow;  that  is 

NQ 

’£A(j,k)a(k)=-ni(j)Vr  (2.28) 

k= 1 

where  NQ  is  the  number  of  elements  used  to  approximate  the  body  surface.  EQ.  2.28 
must  be  satisfied  at  each  null  point.  Repeating  EQ.  2.28  for  each  element  results  in 
the  set  of  linear  algebraic  equations  that  replaces  the  integral  EQ.  2.13  as 

(A)  {o-}  =  (2.29) 

The  selection  of  null  points  ensures  that  A(j,  k)  =  2n  as  j  =  k  (all  the  diagonal  terms 
of  matrix  A).  The  solution  of  EQ.  2.29  gives  the  source  density  values  r  for  each 
quadrilateral.  Once  these  source  densities  are  obtained,  the  induced  velocities  at  any 
point  away  from  the  body  surface  are  given  by 

NQ 

V,=Z  (2.30) 

i= l 
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The  calculation  of  induced  velocity  V{(j)  due  to  each  element  follows  EQ.  2.17  to 
EQ.  2.19  and  transformed  to  the  global  coordinates  as  was  done  for  the  null  points. 

2.4  Potential  Flow  and  Strain  Paths 

For  penetration  problems,  it  is  convenient  to  use  a  Eulerian  coordinate  system  where 
a  point  on  the  penetrometer  is  chosen  as  the  origin.  The  advantage  of  this  system  is 
that  all  the  field  variables  are  fixed  in  space  with  reference  to  this  origin.  Results  of 
solving  the  potential  flow  problem  as  in  EQ.  2.30  generate  the  induced  flow  velocity 
at  a  point  in  space.  Considering  a  short  time  period  dt,  the  amount  of  displacement 
at  that  point  is 

Ui  =  J  Vidt  ~  Vidt  (2.31) 

during  this  time  period.  EQ.  2.31  implies  that  the  induced  velocities  remain  con¬ 
stant  during  each  time  increment  dt.  In  order  for  this  assumption  to  be  valid,  it  is 
necessary  to  use  a  sufficiently  small  dt  for  each  increment  which  results  in  very  small 
displacements.  Therefore,  even  the  final  accumulated  strains  could  be  large,  the  small 
strain-displacement  relationship  is  sufficiently  accurate  for  each  increment  where 


^ = \<fi* + ty 

(2.32) 

Or  according  to  EQ.  2.31 

i<j  =  +  V^it) 

(2.33) 

where 

u  -  2* 

U'J  dx, 

v  = 

dXj 
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In  local  coordinates,  the  induced  strain  tensor  due  to  each  quadrilateral  is 

iij  =  +  Ijjdt )  (2.34) 

where 


/ 

dxj 

The  derivatives  are  determined  analytically  from  EQ.  2.17  to  EQ.  2.19.  Transforming 
to  global  coordinates 

£ij  =  (2.35) 


Summing  the  components  from  all  the  elements,  the  induced  strain  tensor  increment 
becomes 


NQ 

£ij  =  E  e<j(*) 


(2.36) 


fc=i 


Integration  of  EQ.  2.31  is  analytically  intractable.  This  is  basically  because  the  po¬ 
tential  flow  analysis  computes  flow  velocities  only.  To  keep  track  of  the  displacements 
and  strains  in  soil  elements,  at  a  given  time  t,  they  are  integrated  numerically  by  fol¬ 
lowing  a  stream  line  from  a  starting  point  far  ahead  of  the  penetrometer  (Levadoux 
and  Baligh,  1980). 


A ’T  NT 

Dx  =  £  U(m)  =  Vi(m)dt(m)  (2.37) 

m=l  m=l 

NT 

ea  =  E  Cii(m)  (2.38) 

m= 1 

where 

NT  =  number  of  time  increments 
dt  =  duration  of  each  time  increment 
D,  =  displacement  components 
£ij  =  accumulated  strain  tensor 
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The  accumulated  strain  from  EQ.  2.38  could  be  very  large  as  will  be  demonstrated 
later.  However,  it  should  not  be  considered  the  same  as  that  resulted  from  computa¬ 
tions  using  large  strain  theories  (Chen  and  Saleeb,  1982). 


2.5  Two-Dimensional  Panel  Methods 

The  implementation  of  the  above  three-dimensional  panel  method  is  very  time  con¬ 
suming  and  demands  large  computer  memories.  This  is  especially  true  when  a  large 
number  of  panels  are  involved.  The  computation  time  can  be  reduced  by  one  order  of 
magnitude,  if  the  problem  is  simplified  into  2-D.  For  example,  if  a  flat-plate  penetrom¬ 
eter  is  infinitely  wide  in  £3  direction,  then  all  velocity  and  strain  components  related 
to  x3  coordinate  are  zero.  The  2-D  simplification  is  accomplished  by  treating  the 
panels  as  infinitely  long  £3.  In  cases  of  axisymmetric  bodies,  two-dimensionalpanels 
are  a  series  of  frustums  of  cones.  The  induced  velocities  are  derived  by  modifying  the 
integral  EQ.  2.16.  These  procedures  are  well  established  in  aeronautical  engineering 
(Hess  and  Smith,  1962)  and  can  be  easily  formulated  to  compute  strain  paths.  The 
two-dimensionalSPM  programs  have  also  developed  in  this  research  which  are  capable 
of  analyzing  strain  paths  for  infinitely  wide  and  axisymmetric  penetrometers. 


2.6  Evaluation  of  the  panel  method 

A  series  of  preliminary  strain  path  analyses  were  performed  for  a  cone  and  a  flat 
dilatometer  to  validate  the  new  SPM  technique.  Figure  2.4  shows  the  distribution 
of  panels  on  the  cone  and  dilatometer  tips.  The  simulated  cone  was  820  mm  long, 
had  a  60  degree  apex  and  a  diameter  of  35.6  mm.  The  flat  dilatometer  (DMT) 
was  simulated  as  14  mm  thick,  94  mm  wide  and  235  mm  long.  The  penetration 
rate  was  20  mm/sec  in  both  cases.  A  reference  radius  R  is  used  in  presenting  the 
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data.  The  radius  of  cone  (17.8  mm)  and  half  thickness  (7  mm)  of  the  DMT  was 
used  as  the  R  value  for  the  respective  penetrometer.  For  the  cases  to  be  presented, 
the  integration  of  displacements  and  strains  started  at  an  initial  coordinate  of  (0,- 
200R,R)  which  was  200  R  ahead  of  the  penetrometer  tip  (see  Figure  2.4)  for  the 
coordinate  system).  The  initial  integration  interval  was  5R  and  reduced  as  the  to 
penetrometer  distance  shortened.  A  minimum  integration  interval  of  0.02R  was  used 
and  maintained  when  a  distance  of  10R  was  reached.  These  parameters  are  similar 
to  those  used  by  Levadoux  and  Baligh  (1980).  Further  reduction  of  these  integration 
intervals  showed  no  significant  change  in  results. 

In  presenting  the  strain  paths,  the  is, -space  (Prevost,  1978)  is  used,  where 


t\  —  £22 

(2.39) 

£2  =  —£=(£11  —  £33) 

(2.40) 

£3  =  2, 

(2.41) 

04 

u7 

11 

(2.42) 

£s  -  7f’3 

(2.43) 

Compressive  strain  is  negative.  The  Et  space  provides  a  simple  method  to  present  a  six 
component  strain  tensor  in  a  three-dimensional  space.  The  components  are  directly 
related  to  some  major  shearing  modes.  The  £i-strain  is  produced  by  a  shearing  mode 
of  triaxial  tests,  the  ^-strain  is  caused  by  lateral  expansion  and  E3,  E4  and  F^-strain 
is  produced  by  direct  simple  shear  tests. 

A  total  of  576  panels  were  used  to  cover  the  flat  DMT  which  has  a  rather  simple  geom¬ 
etry.  The  adequacy  of  this  panel  configuration  was  first  evaluated  in  two-dimensional. 
The  DMT  was  considered  infinitely  long  in  aq  direction  shown  in  Figure  2.4  .  Panels 
behind  the  DMT  tip  are  parallel  to  the  direction  of  penetration  and  their  effects  to 
strain  paths  are  insignificant.  Figure  2.5  shows  the  strain  paths  in  terms  of  Ei  versus 
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(i)  Cut 


Figure  2.5:  Effects  of  resolution  on  strain  paths. 


Figure  2.6:  Induiced  strain  paths. 
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Displacement/R  in  x%  direction 


Figure  2.7:  Induced  displacements. 
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E3  resulted  from  three  analyses  using  different  number  of  elements.  The  first  curve 
(32  elements)  had  half  the  number  of  panels  as  that  of  42  elements  on  the  DMT  tip. 
The  last  curve  (60  elements)  had  double  the  number  of  panels  on  the  tip.  We  can  say 
that  these  three  curves  are  identical,  considering  the  qualitative  nature  of  the  strain 
path  method,  even  as  the  resolution  of  panels  changed  by  four  fold.  Further  increase 
of  the  panel  number  resulted  in  no  significant  differences  in  strain  paths.  Because 
the  DMT  has  a  uniform  shape  in  the  ij  direction,  the  number  of  panels  in  this  di¬ 
mension  has  little  effect  to  the  results.  The  configuration  shown  in  Figure  2.4  was 
selected  basically  to  give  a  reasonable  length  to  width  ratio  for  each  panel.  A  total 
of  2880  elements  were  used  to  simulate  the  cone  (Figure  2.4)  in  a  three-dimensional 
SPM  computation.  This  is  far  more  than  necessary  to  give  reasonable  results.  The 
higher  element  number  was  used  primarily  to  demonstrate  that  the  round  surface  of 
a  cone  can  be  realistically  represented  by  a  large  number  of  such  plane  quadrilater¬ 
als.  It  should  also  be  noticed  that  there  was  no  need  to  avoid  the  sharp  corner  at 
the  base  of  the  cone  tip  as  in  other  methods  (Baligh,  1985)  of  strain  path  compu¬ 
tations.  The  analysis  was  also  performed  using  a  technique  based  on  axisymmetric, 
two-dimensionalpanel  method  where  the  cone  surface  was  represented  by  25  line  ele¬ 
ments.  The  displacement  and  strain  components  (Figure  2.6  and  Figure  2.7)  of  the 
cone  from  both  the  two-dimensionaland  three-dimensional  analyses  are  essentially 
identical  to  those  reported  by  Levadoux  and  Baligh  (1980). 
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Chapter  3 

The  Strain  Field 


The  new  numerical  technique  enabled  strain  path  analyses  be  performed  for  flat- 
plates  which  were  not  possible  before.  Analyses  were  conducted  for  the  DMT,  K0- 
Stepped  Blade  and  continuously  tapered  blade  for  the  first  time,  using  the  new  three- 
dimensional  strain  path  technique.  Dimensions  of  these  penetrometers  used  in  the 
computations  are  shown  in  Figure  3.1.  This  chapter  describes  the  characteristics  of 
the  strain  field  during  a  flat-plate  penetration. 

A  slightly  different  coordinate  system  was  adopted  herein  after.  Figure  3.2  shows  the 
coordinate  system  used  and  its  relation  to  the  flat-plate.  All  linear  dimensions  are 
normalized  with  respect  to  a  reference  radius  R.  For  axisymmetric  penetrometers  R 
is  the  actual  radius,  whereas  R  equals  the  half  thickness  of  the  blade  for  flat-plates 
Figure  3.1.  Computations  of  soil  displacements  and  strains  follow  a  stream  line  started 
far  ahead  of  the  penetrometer,  typically  at  y  =  —200 R  as  described  previously.  In 
presenting  the  strains,  the  starting  coordinates  (x,z)  of  the  associated  stream  line  are 
given  and  the  F,  space  (Prevost,  1978)  as  mentioned  before  is  used.  A  strain  ratio 
SRi  is  used  to  indicate  the  relative  importance  for  each  strain  component  during 
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Figure  3.1:  Dimensions  of  the  flat-plates  analyzed. 


penetration. 


5/2,  = 


EL,  £,2 


Plots  of  E ,  and  5/2,  versus  y/R  describe  how  they  vary  as  the  soil  moves  relative  to 
the  penetrometer. 


Strain  paths  ( E2  versus  £3)  and  soil  displacements  following  streamlines  started  at 
z  =  1/2,  3R,  5R  and  10R,  and  x  =  OR  are  shown  in  Figure  3.3.  The  results  indicate 
that  the  magnitudes  of  strains  and  displacements  decrease  with  distance  to  the  plate. 
However,  the  pattern  of  strain  paths  as  shown  in  Figure  3.3a  remains  similar.  As  a 


reference,  the  strain  path  during  a  cone  penetration  (CPT),  along  a  stream  line  started 
at  z  =  1/2  and  x  =  OR  is  also  included  in  Figure  3.3.  Both  penetrometers  cause 
significant  strain  reversals  as  demonstrated  in  Figure  3.3a.  However,  an  important 
difference  between  the  CPT  and  the  DMT  is  the  much  higher  magnitude  of  /^-strain 
which  occurs  during  the  CPT  penetration.  This  is  the  result  of  a  very  large  tensile 
circumferential  strain  which  develops  during  the  CPT  penetration.  Plots  of  strain 
ratios  following  stream  lines  at  z  =  1/2  and  x  =  OR  (Figures  3.4  and  3.5)  show  that 
in  both  cases,  E\  is  the  predominant  strain  component  before  the  soil  particle  reaches 
the  tip  of  penetrometers  (y  =  0/2).  For  the  DMT,  £3  and  E2  exchange  their  role 
as  the  predominant  strain  component  after  passing  the  tip  with  £2  becoming  the 
predominant  component  when  approaching  the  flexible  membrane  (13/2  <  y  <  17/2). 
Beyond  the  tip  (as  y  >  0/2),  E2  becomes  the  major  strain  component  in  the  case  of 
CPT  penetration,  which  provides  some  justification  to  use  cylindrical  cavity  expansion 
theory  in  the  interpretation  of  CPT  results. 


3.1  Three  Dimensional  Effects 

To  demonstrate  the  three-dimensional  nature  of  the  DMT  penetration,  Figures  3.6 
and  3.7  show  a  quadrant  of  the  contours  of  displacements,  ds  and  shear  octahedral 
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strains,  7 oct  at  y  =  13 R.  This  is  approximately  the  same  level  at  the  center  of  the 
DMT  flexible  membrane.  The  displacement  is  computed  as: 


ds  = 


(. dx 2  +  dz2)xl 2 
R 


and  the  shear  octahedral  strain 


(3.2) 


2  JL 

7 oct  =  g  [(£xx  ~  £yy)  T  (Sxx  —  £zz)  T  {&yy  ~  &zz)  "F  fl(^xv  £yz  T  ^zx)]  (3-3) 

The  magnitudes  of  shear  octahedral  strains  are  much  higher  towards  the  corner  of  the 
plate  (Figures  3.6  and  3.7).  The  limited  plate  width  affects  not  only  the  magnitude 
but  also  the  characteristics  of  the  strains.  Figure  3.8  shows  £2  versus  £3  following 
stream  lines  at  x  =  OR,  3.4 R  and  6.7 R  ( z  —  1 R  for  all).  At  x  =  6.7 R,  the  stream 
line  passes  the  corner  of  the  DMT.  Along  the  center  of  the  DMT  (i.e.  at  x  =  OR) 
£XI  is  compressive  while  £zz  is  tensile.  The  £2-strain  increases  towards  the  corner  of 
the  DMT  and  approaches  the  magnitude  which  occurs  during  CPT  penetration  or  a 
cylindrical  cavity  expansion  Figure  3.3a).  Immediately  adjacent  to  the  DMT  corner, 
there  is  a  reversal  of  signs  for  both  exx  and  ezz.  This  results  in  a  positive  ^-strain 
along  the  corner  of  the  blade  as  shown  in  Figure  3.8. 


To  further  demonstrate  the  three-dimensional  nature  of  the  DMT  penetration,  the 
corresponding  shear  octahedral  strain  and  displacement  values  from  2-D  (w/t  =  00) 
analyses  are  included  in  Figures  3.6  and  3.7.  The  results  clearly  indicate  that  the  2-D 
simplification  would  significantly  underestimate  the  shear  octahedral  strains.  Also, 
the  strain  and  displacement  components  attenuate  much  slower  with  distance  in  2-D 
analyses. 


A  series  of  strain  path  analyses  were  performed  using  different  widths  for  the  blade  to 
investigate  the  significance  of  three-dimensional  effects.  Figure  3.9  shows  the  strain 
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paths  in  terms  of  E 2  and  ^-strains  following  stream  lines  at  x  =  OR  and  z  =  \R. 
In  these  analyses,  the  width  over  thickness  ratios,  w/t  varied  from  2  to  infinity. 
For  a  blade  with  infinite  width,  the  analysis  actually  produces  a  2-D  solution.  The 
results  presented  in  Figure  3.9  show  that  with  w/t  =  2,  the  .E^-strain  approaches  the 
corresponding  values  obtained  during  CPT  (Figure  3.3a).  As  w/t  increases,  the  strain 
paths  approach  the  2-D  solution.  Above  a  w/t  of  approximately  20,  the  computed 
strain  paths  resemble  that  from  2-D  analysis  where  the  plate  width  was  considered 
infinite.  Further  increase  of  w/t  ratios  beyond  20  did  not  result  in  any  significant 
differences  in  strain  paths.  The  results  of  Figure  3.9  suggest  that  the  strain  paths 
immediately  adjacent  to  the  center  (x  =  OR)  of  a  penetrating,  thin  flat-plate  would 
approximate  the  2-D  case. 


3.2  Effects  of  Tip  Apex  Angles 

Using  the  2-D  analyses,  it  was  possible  to  evaluate  the  effects  of  the  DMT  apex  angles 
without  the  complication  caused  by  the  limited  plate  width.  Figure  3.10  shows  the 
predicted  soil  distortions  around  the  penetrometer  with  four  different  tip  apex  angles 
ranging  from  10  to  60  degrees.  The  DMT  has  an  equivalent  apex  angle  of  approxi¬ 
mately  20  degrees.  The  strain  paths  along  stream  lines  at  z  =  IR  show  that  £3-strain 
is  very  sensitive  to  the  apex  angle  Figure  3.11.  At  an  apex  angle  of  60  degrees,  £3 
is  much  larger  than  £2  which  is  not  the  case  during  a  comparable  CPT  penetration 
(see  Figure  3.3a).  The  predominant  strain  component  during  a  flat-plate  penetration 
is  therefore,  affected  by  the  tip  apex  angle.  Further  results  to  be  shown  later  would 
demonstrate  that  the  same  is  true  from  actual  three-dimensional  strain  path  analyses. 

The  contrast  in  soil  distortions  between  the  60  degree  plate  and  the  DMT  are  com¬ 
parable  to  that  presented  by  Levadoux  and  Baligh  (1980)  for  cones  with  60  and  18 
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degree  apex  angles.  One  might  conclude  that  if  the  DMT  penetration  does  induce 
less  soil  disturbance  than  the  CPT  as  stated  by  Marchetti  (1981),  it  is  most  likely  due 
to  the  sharper  apex  angle  of  the  DMT  and  not  necessarily  because  of  the  geometry 
(i.  e.  ,  flat  versus  cylindrical). 


3.3  Effects  of  Geometry 

From  the  context  of  the  interpretation  theories,  Handy,  et  al.  (1982)  actually  assumed 
that  E2  is  the  predominant  strain  component  during  penetration  and  E2  is  propor¬ 
tional  to  the  blade  thickness.  The  Stepped-Blade  was  made  thin  enough  so  that  the 
penetration  can  be  considered  plane  strain,  which  is  a  requirement  for  Eq.  ??  to 
be  valid.  In  an  attempt  to  increase  the  number  of  stress  measurements  and  improve 
the  reliability  of  stress  extrapolation  a  continuously  tapered  blade  was  conceived  as 
a  taxadjunct  to  the  Stepped  Blade  (Handy  and  Lutenegger,  1985).  The  development 
of  this  device  by  the  second  author  has  yielded  a  blade  which  has  a  continuous  taper 
of  3  degrees  and  is  instrumented  with  seven  stress  sensors  along  one  face.  Currently, 
the  interpretation  of  results  from  the  tapered  blade  has  followed  the  extrapolation 
procedure  proposed  by  Handy,  et  al.  (1982). 

Figures  3.12  through  3.16  show  the  results  of  the  three-dimensional  analyses  for  the 
/f0-Stepped  and  tapered  blade  penetrometer  in  terms  of  displacements,  E,  strains  and 
strain  ratios  SR*.  Ail  computations  follow  stream  lines  started  at  x  =  OR,  y  =  —200 R 
and  z  =  IR.  These  strain  paths  are  considered  representative  in  the  area  adjacent  to 
the  pressure  sensors  (diameter  less  than  25.4  mm)  on  the  blades. 

For  the  /f0-Stepped  Blade,  the  results  show  that  displacements  in  the  z  direction 
are  approximately  proportional  to  the  blade  thickness  (Figure  3.12).  However,  the 
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strain  components  are  far  from  those  assumed  by  Handy,  et  al.  during  a  stepped- 
blade  penetration.  According  to  Figures  3.12  and  3.13,  £3-strain  is  the  predominant 
strain  component  and  increases  with  the  blade  thickness.  The  i?2-strain  remains  ap¬ 
proximately  constant  and  appears  independent  of  blade  thickness.  Also,  all  strain 
components  undergo  significant  strain  reversals  between  successive  steps.  The  strain 
paths  imply  that  the  stresses  measured  on  the  face  of  a  stepped-blade  are  equiva¬ 
lent  to  the  normal  stress  response  during  a  cyclic  simple  shear  test  where  2i3-strain 
is  the  major  component.  Since  the  width  is  constant,  the  AVStepped  Blade  has 
a  w/t  ratio  which  starts  at  about  22  at  the  first  step  and  decreases  to  8.7  at  the 
fourth  step  (Figure  3.1).  In  view  of  the  aforementioned  w/t  studies  for  the  DMT 
(Figure  3.9),  the  stepped-blade  penetration  is  not  likely  to  be  plane  strain,  at  least 
for  the  thicker  part  of  the  blade  due  to  the  less  favorable  w/t  ratios.  This  was  in  part 
recognized  by  Handy,  et  al.  (1982)  who  hypothesized  that  thicker  blades  would  re¬ 
sult  in  a  transition  to  cylindrical  cavity  behavior  and  thus  produce  a  "limit”  pressure. 

The  soil  disturbance  induced  by  the  tapered  blade  penetration  (Figures  3.12,  3.15  and 
3.16)  is  less  erratic  compared  with  that  of  the  stepped-blade.  As  in  stepped-blade, 
the  £3-strain  is  the  predominant  component.  After  passing  the  tip,  all  three  strain 
components  remain  relatively  constant  as  the  blade  thickness  increased  from  6.6  mm 
to  12  mm.  There  are  no  strain  reversals  as  there  are  no  "steps”  and  thus  the  soil  does 
not  undergo  an  unloading  sequence. 
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Figure  3.2:  The  coordinate  system. 
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placement/!?  in  y— direction 


0.0  0.2  0.4  0.6  0.8  1.0 


Displacement/!?  in  z-direction 

Figure  3.3:  Strain  paths  and  displacements  during  DMT  penetration. 
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Figure  3.4:  Strain  ratios  during  CPT  penetrations. 


Figure  3.7:  Contours  of  shear  octahedral  strains. 


Figure  3.8:  Strain  paths  approaching  the  corner  of  DMT. 
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Figure  3.9:  Effects  of  w/t  on  strain  paths. 
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(b)  DMT 


(d)  60* 


Figure  3.10:  Effects  of  tip  apex  angle  on  soil  disturbance. 
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Figure  3.11:  Effects  of  tip  apex  angle  on  strain  paths 
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Displacement /R  in  z- direction 
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Figure  3.14:  Strain  ratios  during  the  stepped  blade  penetration. 
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Figure  3.16:  Strain  ratios  during  the  tapered  blade  penetration. 


Chapter  4 

Deviatoric  Stress  Model 


The  strains  are  completely  decoupled  from  stresses  in  the  strain  path  computations. 
Therefore,  the  strain  path  analysis  itself  does  not  yield  any  information  about  stresses. 
In  an  attempt  to  improve  the  situation,  Baligh  (1985)  proposed  to  relate  stresses  to 
strains  through  an  inverted  stress-strain  model.  Most  of  the  conventional  stress- 
strain  models  are  so  called  ” stress-space"  models  where  stresses  are  the  independent 
variables.  When  inverted,  there  is  no  direct  mechanism  to  relate  deviatoric  strains 
(volumetric  strains  =  0  during  undrained  penetration)  to  normal  stresses.  Thus,  only 
the  deviator  stresses,  s,j  can  be  estimated  as 

sij  =  /(£«>)  (M) 

Ideally,  the  soil  model  used  in  estimating  the  stress  field  should  be  able  to  consider 
a)  complicated  strain  paths  including  strain  reversals  and  large  plastic  strains;  b) 
initial  and  stress  induced  anisotropy;  c)  post  peak  behavior;  and  d)  time  dependent 
deformations  (e.g.,  undrained  crep  relaxation).  Such  a  compreher  sive  model  is  not 
available  at  the  present  time.  Many  existing  models  can  deal  with  complicated  strian 
paths  and  anisotropy.  Few  models,  however,  fares  well  in  handling  post  peak  behavior 
and  time  dependent  deformation.  A  strain  based  Modefied  Cam  Clay  (MCC)  Model 
was  selected  to  estimate  the  corresponding  stresses  from  the  strains  generated  in  the 
SPM  analyses.  This  selection  was  based  on  the  following  reasons: 
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1.  MCC  is  an  isotropic  model  capable  of  handling  the  complicated  strain  paths 
as  demonstrated  by  Kavvadas(19S2).  Because  of  the  large  strains  involved,  a 
sophisticated  model  is  likely  to  mask  out  important  characteristics  predicted 
by  SPM.  In  addition,  the  SPM  generated  strain  field  should  be  considered  as 
qualitative  because  of  the  way  it  is  computed.  It  is  hardly  justifiable  to  consider 
issues  such  as  anisotropy  and  time  dependent  deformations; 

2.  MCC  model  is  simple  in  establishing  the  input  parameters  and  it  has  been  used 
extensively; 

3.  Contrast  to  cone  penetrations,  it  is  not  known  if  shear  stress  induced  pore  water 
pressure  is  significant  in  comparison  wi  h  the  mean  normal  stress  related  pore 
water  pressure.  Because  MCC  is  an  effective  stress  model,  the  shear  stress 
induced  pore  water  pressure  is  computed  directly. 

The  model  was  calibrated  according  to  a  series  of  controlled  gradient  consolidation 
and  consolidated  undrained  triaxial  (CIU)  tests.  This  chapter  describes  the  imple¬ 
mentation  of  the  MCC  model  using  strains  as  the  independedent  variables  following 
the  development  by  Kavvadas  (1982). 


4.1  General  Concepts  of  Plasticity 


It  is  common  in  conventional  plasticity  to  divide  the  total  strain  increment,  corre¬ 
sponding  to  an  infinitesimal  stress  increment,  into  two  parts: 


1.  the  reversible  portion  called  the  elastic  strain  increment,  £'•,  and; 

2.  the  reversible  plastic  strain  increment,  i?- 


Thus, 


£u  — 


£  +  £y 
•J  ~  u 


(4.2) 


The  conventional  definition  of  deviatoric  and  volumetric  parts  of  the  strain  increment 
are  accepted.  They  can  be  stated  as: 


£  —  £ijl  ij 


(4.3) 
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and, 


(4.4) 


4.1.1  Elastic  Part  of  Elastoplastic  Deformation 

For  most  loading  sequences  the  elastic  part  of  soil  deformation  is  much  smaller  than 
the  plastic  part.  Thus,  simple  isotropic  (linear  or  nonlinear)  models  are  justified  and 
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are  commonly  used  (e.g.,  Roscoe  and  Burland,  1968;  Davis  and  Mullenger.  1978). 
The  following  isotropic  relationship  with  two  independent  parameters  is  used. 

W  =  Kie 

i,i  =  2  Ge^  (4.10) 

where,  K  is  the  elastic  bulk  modulus  and  G  is  the  elastic  shear  modulus. 


For  normally  consolidated  clays  that  not  sensitive: 

,,  1  +  eo_  , .  . , , 

A  = - a  (4.1 1 ) 

K 

where  e0  is  the  initial  void  ratio  and  k  is  a  materia]  constant  and  equals  the  slope  of 
swelling  line  in  the  Cam-clay  model. 


The  elastic  shear  modulus  G  affects  the  initial  stiffness  and  also  the  stiffness  of  the 
unloading  part  of  the  shear  stress  strain  curve.  The  Modified  Cam-Clay  model  ne¬ 
glects  elastic  shear  strain  which  means  G  is  infinitely  large.  Such  an  assumption  is 
not  realistic  and  leads  to  numerical  problems.  Therefore.  Kavvadas  (1982),  proposed 
that  the  ratio  |  be  a  constant.  The  value  of  this  constant  is  dictated  by  elastic 


Poisson’s  ratio,  v. 

2  G  „  1-217 
K  1  -+■  17 

Thus  I\  and  G  are  both  functions  of  the  current  level  of  a. 
elastic  portion  of  soil  deformation  is  non-linear  (and  isotropic). 


(4.12) 


This  means  that  the 


4.1.2  Plastic  Part  of  Elastoplastic  Deformation  of  Soil 

The  relationship  between  the  plastic  part  of  strain  increment  and  effective  stress 
increment  can  be  stated  in  a  genera!  form  as: 

iPH  =  APU  (4.13) 

where,  A  is  a  scalar  controlling  the  the  magnitude  of  strain  increment  and  is  a 
second  order  tensor  that  determines  the  direction  of  strain  increment.  Further,  PX} 
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is  independent  of  a.  A  can  be  written  as  a  linear  combination  of  the  components  of 
effective  stress  increment. 

A  =  £(<?, A)  (4.14) 

where,  H  is  the  elastoplastic  modulus  and  QtJ  is  a  second  order  tensor  that  deter¬ 
mines  the  magnitude  of  incremental  strain  through  A.  The  second  order  tensors  Pt] 
and  Qij,  like  the  strain  tensor  can  be  divided  into  their  volumetric  components  (P, 
Q)  and  deviatoric  components  (P-j,  Q[j)- 

Equation  4.13  can  be  expressed  separately  for  volumetric  and  deviatoric  components 
in  the  form: 


ep  =  A  P 


=  Api) 


(4. 15) 


combining  equations  4.14  and  4.15  wre  get: 

A  =  +  Q&)  (4.16) 

Finally  equations  4.13  and  4.14  lead  to: 

efj  =  jjiQkiVkriPij  (4.17) 

Soil  penetration  problems  such  as  cone  penetration  and  Dilatometer  penetration  are, 
inherently,  strain  controlled.  The  strains  induced  around  the  penetrometer  are  a  re¬ 
sult  of  kinematic  constraints.  Thus,  the  problem  is  posed  inversely:  Given  the  applied 
strain  increment  it  is  required  to  establish  the  corresponding  stress  field  around  the 
penetrometer.  The  stress  field  comprises  of  the  effective  stress  increment  atJ  or  its 
component  a  and  s,j. 


Combining  equations  4.2  and  4.10  we  get: 


t  =  Kie  =  K(i  -  ep) 
■  2 Gi\i  =  2G(ei}  -  cf,) 


(4. IS) 
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using  equations  4.10  and  4.15 


ir  =  I<(i  -  A P) 
sn  =  2 GUij  -  A P(j) 

Now,  using  equations  4.15  and  4.19 


(4.19) 


//A  =  +  QtJ1f  =  2 GQUeij  -  AP')  +  KQ(i  -  A P) 


(4.20) 


Solving  for  A; 


A  = 


2G(Q;.,eIJ)  +  KQe 


H  +  2G(Q[JP(j)  +  KQP 


(4.21) 


Equation  4.18  gives  the  increment  of  mean  effective  stress  and  also  the  deviatoric 
stress  increment  in  terms  of  volumetric  and  deviatoric  strain  increments.  Equa¬ 
tion  4.21  is  linear  in  stress  and  strain  increments.  A  is  a  linear  function  of  strain 
increments,  therefore  if  A',  G ,  P,  P(j,  Q,  Q\-  and  H  are  known  the  above  equation 
can  be  implemented  in  a  numerical  scheme. 


4.2  The  Yield  Function  and  Flow  Rule 

The  MCC  model  assumes  an  ellipsoidal  yield  surface,  in  the  effective  stress  space, 
symmetric  about  its  hydrostatic  axis  (Figure  4.1a).  The  effective  stress  space  is 
comprised  of  the  deviatoric  stress  subspace  stJ  and  the  mean  effective  stress  a.  The 
yield  function  for  such  a  yield  surface  is  written  as; 

f  (<j ,  a,  Sij)  SijS{j  c  cr(2a  <r) 

or, 

f{a,a,s,})  =  r2  -  c2a{ 2a  -  a)  (4.22) 

where, 

•  r  =  {sijSij)i  =  V2J2 
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Figure  4.1:  Yield  surface  of  MCC  model  in  stJ  —  a  space. 
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•  J2  is  the  second  invariant  of  deviatoric  stress  tensor 

•  c  is  the  material  constant  and  is  equal  to  the  ratio  of  the  axes  of  the  yield  surface 

*=\Hm 

•  M  =  q/p  and,  q  =  (cri  —  07/),  p  —  a  at  critical  state. 

•  a  is  a  state  variable  called  the  hardening  parameter.  It  controls  the  size  of  the 
yield  surface.  Changes  of  a  during  deformation  allow  the  strength  to  increase 
with  consolidation. 

4.2.1  Plastic  Loading  and  Elastic  Unloading  from  Elasto- 
plastic  States 

The  flat-plate  penetration  mechanism  involves  complex  deformation  in  the  soil  el¬ 
ement  under  consideration.  Certain  stress  components  decrease  while  other  stress 
components  increase,  therefore  the  definitions  of  loading  and  unloading  have  to  be 
unambiguous.  The  MCC  model  satates  that: 

the  application  of  an  infinitesimal  effective  stress  increment  represents  plastic  loading 
if  both  plastic  and  elastic  strains  occur,  and  represents  elastic  unloading  if  only  elastic 
strains  occur. 

Therefore  starting  from  a  stress  state  on  the  yield  surface  i.e.  normally  consolidated 
state,  elastic  unloading  implies  a  retreat  of  the  stress  point  from  the  yield  surface 
causing  zero  plastic  strain.  Whereas,  plastic  loading  implies  an  effective  stress  incre¬ 
ment  directed  outside  the  yield  surface  thus  leading  to  a  new  elastoplastic  state.  A 
smooth  transition  between  the  elastic  unloading  and  plastic  loading  at  the  boundary 
is  effected,  in  MCC  model,  by  the  limiting  case  assumption  that;  an  effective  stress 
increment  at  a  tangent  to  the  yield  surface  causes  no  plastic  strains.  This  requirement 
is  satisfied  if  the  second  order  tensor  Q,:  is  defined  as  the  gradient  of  the  yield  surface 
at  the  stress  point  representing  the  current  state  of  stress. 


Decomposing  Qij  into  its  volumetric  and  deviatoric  parts  we  get: 


(4.24) 


il  _i,M. ,  u 

dsh  3  {dst,Iu]!ii 


(4.25) 


Having  defined  Qij  as  the  gradient  of  the  yield  surface  at  the  stress  point  representing 
the  current  state  of  stress,  the  MCC  model  proposes  a  quantitative  relation  to  distin¬ 
guish  between  the  plastic  loading  and  elastic  unloading  from  elastoplastic  states.  For 
the  case  of  elastic  unloading  the  plastic  strain  increment  is  defined  to  be  zero  and  the 
stress  point  retreats  from  the  yield  surface.  Since  Qij  is  directed  towards  the  outer 
normal  to  the  yield  surface  it  follows  that  for  elastic  loading: 


QijVij  <  0  (4.26) 

using  the  decomposition  of  the  effective  stress  increment  and  Qij  into  its  volumetric 
and  deviatoric  parts  we  can  show  that: 

Q'ijSij  +  Q$<  0  (4.27) 

Using  equations  4.10  and  4.27  the  following  criterion  for  elastic  unloading  is  obtained: 

2 G{Q’ijtij)  +  KQe  <  0  (4.28) 

and  the  conjugate  criterion  for  plastic  loading: 

2G{Q'ijix])  + KQi  >  0  (4.29) 


Having  defined  the  plastic  loading  and  elastic  unloading  criteria  the  components  of 
Qij  can  be  specified  further  using  equations  4.23,  4.24  and  4.25  for  the  specific  yield 
function  of  the  MCC  model: 


<5  =  2c2(cr  -  or) 


2s 


u 


(4.30) 

(4.31) 
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Furthermore  the  model  employs  an  associated  flow  rule,  so  that: 


P  =  Q 

PI,  =  <?;,  (4.32) 

This  means  that  the  direction  of  the  plastic  strain  increment  (given  by  )  coincides 
with  the  normal  to  the  yield  surface  (given  by  QtJ). 


4.3  Hardening  Characteristics 

In  general,  the  hardening  characteristics  of  a  model  describe  the  evolution  of  model 
parameters  specifying  the  shape,  size  and  orientation  of  the  yield  surface  during  plas¬ 
tic  deformation.  The  MCC  model  utilized  herein  adopts  an  isotropic  hardening  rule 
and  a  is  the  only  parameter  describing  the  yield  surface  of  the  MCC  model.  The 
hardening  rule  determines  the  evolution  of  a  during  plastic  deformation  of  soil.  This 
allows  the  yield  surface  to  expand  in  the  stress  space  without  changing  its  shape. 


The  MCC  model  assumes  that  for  one-dimensional  consolidation  or  isotropic  con¬ 
solidation  of  normally  consolidated  clay  the  stress  path  gives  a  linear  relationship 
between  the  void  ratio  e  and  the  natural  logarithm  of  mean  effective  stress  a,  with 
slope  A: 

A=  de 
d(ln(T ) 

A  stress  path  for  the  above  mentioned  one-dimensional  consolidation  or  isotropic 
consolidation  of  normally  consolidated  clay  is  also  called  radial  stress  path.  The 
MCC  model  assumes  that  for  plastic  loading  along  a  radial  stress  path  (in  effective 
stress  space)  the  hardening  parameter  a  is  proportional  to  to  the  mean  effective  stress 
a  and  thus, 


dev  =  de  —  dee  =  —X.d(ln(r)  +  Kd(lno) 
dep  =  —(A  —  rc)d(lna)  =  —(A  —  K)d(lna)  (4.34) 
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or: 


ep  =  —(A  —  «)— 


(4.35) 


where,  ep  is  the  plastic  increment  of  void  ratio  and  a  is  the  increment  of  the  size 
parameter  a.  However,  from  the  definition  of  void  ratio: 


ep  =  -(1  +  e0)ip 


(4.36) 


where,  eo  is  the  initial  void  ratio  and  ip  is  the  volumetric  component  of  the  plastic 
strain  increment.  Combining  equations  4.35  and  4.36  we  get: 


1  +  eo_.„ 
a  =  - aev 

A  —  K 

or,  combining  equation  4.15  with  the  above  result  we  get: 

-i.  1  +  eo_;  n 

a  =  - aAP 

A  —  K 


(4.37) 


(4.38) 


4.4  The  Elastoplastic  Modulus  H 

During  plastic  loading  the  stress  point  remains  on  the  yield  surface.  Thus  the  effective 
stress  increment  cr.j  and  the  change  in  the  size  of  the  yield  surface  must,  continuously, 
satisfy  the  condition  /  =  0.  This  is  the  well  known  Prager’s  Consistency  condition 
which  specifies  that: 

/(«,*,,)  =  0  (4.39) 

using  the  chain  rule 

Qi&j  +  !=:«  =  0  (4.40) 

using  equations  4.15,  4.37,  4.38  and  solving  for  H 

U  _  1  4"  e°  I  \  _  n  /  A  Al\ 


A  —  k  da 

using  the  definition  of  the  yield  surface  we  can  obtain: 

^  -  -9,*, 

_ —  —  “C  (X 

oa 

and  using  equation  4.30  and  associated  flow  rule: 


(4.41) 


(4.42) 


P  =  2  c2(ct  —  a) 


(4.43) 


60 


4.5  Effective  Stress  Path  for  Undrained  Deforma¬ 
tion  of  Clay 

Figure  4.1a  is  a  geometrical  representation  of  the  ellipsoidal  yield  surface  in  the  effec¬ 
tive  stress  space  (stJ,  cr).  The  point  R  on  the  yield  surface  represents  an  elastoplastic 
state  of  stress.  Figure  4. IB  represents  the  same  yield  surface  in  the  e  —  a  space,  e  be¬ 
ing  the  void  ratio.  The  model  assumes,  that,  for  a  given  elastoplastic  stress  state  the 
point  A  at  the  tip  of  the  yield  surface  (Figure  4.1a)  is  represented  by  the  intersection 
of  the  virgin  isotropic  consolidation  curve  and  the  isotropic  rebound  curve  (Figure 
4.1b).  It  is  evident  from  Figure  4.1a  that: 


cr  o  =  2  a 


(4.44) 


where  a  is  the  value  of  effective  stress  at  the  tip  of  the  yield  surface. 


Let  us  define  tj  as  the  ratio  of  r  to  the  mean  effective  stress  a: 

r 

The  equation  of  the  yield  surface  (EQ.  3.22)  can  be  expressed  as: 


(4.45) 


a 


<7q  c 2  +  7?2 


(4.46) 


This  is  identical  to  the  equation  of  the  yield  surface  given  by  Roscoe  &  Burland 
(1968),  i.e.  it  describes  an  ellipsoid. 

Undrained  deformation  is  given  by  a  straight  line  parallel  to  the  a  axis  in  the  e  —  a 

space.  Furthermore,  assuming  that  the  virgin  isotropic  consolidation  curve  and  the 

critical  state  curve  are  straight  lines  in  the  e  — a  space,  with  slopes  A  and  k  respectively, 

we  have:  _  __ 

Xln(zr)  =  nln(^r)  (4-47) 

cr  o  Oq 


—  —  i  -  a 

a  a  1  * 


<*t  o  o 


(4.48) 
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combining  equations  3.34  and  3.35  we  get: 


(4.49) 


where,  ce  and  cq  are  the  mean  effective  stresses  associated  with  points  E  and  A  re¬ 
spectively  in  Figure  4.1A  and  4. IB.  The  above  equation  gives  the  effective  stress  path 
during  undrained  shear  of  clay  with  a  void  ratio  e  associated  wit  a  specific  value  of 
the  mean  effective  stress,  ae. 


4.6  Implementation  of  the  MCC  Model 

The  MCC  model  was  used  to  study  a  variety  of  fiat-plate  penetrometers  and  their  two 
dimensional  approximations.  Table  4.1  summarizes  the  relations  needed  to  employ 
the  MCC  Model.  The  model  was  implemented  in  the  following  sequence: 

1.  Examine  the  value  of  a  pre-defined  ’flag’  to  check  whether  the  current  state  of 
stress  of  the  soil  element  is  elastic  or  elastoplastic.  If  the  state  is  elastic  go  to 
(5),  else  continue; 

2.  Equations  4.28  and  4.29  is  the  criterion  used  to  distinguish  between  plastic 
loading  and  elastic  unloading.  If  the  strain  increment  causes  elastic  unloading 
go  to  (5),  else  continue; 

3.  Compute  P,j,  H ,  A.  Use  the  elastoplastic  stress  strain  increment  relation  4.18 
to  calculate  the  effective  stress  increment  corresponding  to  the  applied  strain 
increment; 

4.  The  effective  stress  state  is  updated.  Then,  the  hardening  rules  are  employed 
to  update  the  remaining  state  variables  are  updated.  Exit  the  algorithm; 

5.  This  step  in  reached  if  the  state  of  the  soil  element  is  elastic  or  if  elastic  unloading 
takes  place.  The  elastic  incremental  stress-strain  relations  are  used  in  order 
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Incremental  stress-strain  relations  for  the  MCC  model 

k  =  0.434C’, 

c=  y/2/3 M 

A  =  0.434Cc 

K  —  l±sa 

K 

£  —  £{j  Itj 

e,j  =  iij  -  \ilij 

Q  =  P  =  2c2(<t  —  a) 

Q'ij  =  P'a  =  2-0- 

H  =  2c2Wa(^)P 

A  2G(Q,,e,,)+K<3£ 

A  ~  h+2G(<?'jp;j)+/v<3P 

a  =  K  (i  —  A  P) 

in  =  2 G(e«  -  A/V) 

t  =  (^)aAP 

a  =  a  +  <7 

a  =  a  +  a 

$xj  —  $ij  “f-  Sfj 

Table  4.1:  Incremental  stress-strain  relations  for  MCC  model 
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to  compute  the  effective  stress  increment  corresp*  iding  to  the  applied  strain 
increment; 

6.  The  new  effective  stress  state  aij  is  then  examined  to  see  if  it  is  still  inside  the 
yield  surface  or  on  the  yield  surface,  or  if  it  has  crossed  the  yield  surface.  The 
"’•st  two  cases  require  the  appropriate  setting  of  the  ’flag’  variable.  However, 
if  the  stress  state  has  crossed  the  yield  surface  then  it  means  that  the  strain 
increment  is  too  large  and  the  following  step  must  be  followed; 

7.  The  point  of  intersection  is  then  calculated  and  the  fraction  of  the  strain  incre¬ 
ment  responsible  for  this  stress  state  is  computed; 

8.  The  stress  state  is  then  updated  using  the  point  of  intersection  of  the  with  the 
yield  surface  as  the  new  stress  state.  The  remaining  strain  increment  is  then 
evaluated  subtracting  the  strain  increment  required  to  reach  yield  surface  from 
the  applied  strain  increment.  The  flag  is  appropriately  set  and  algorithm  control 
is  transferred  to  (3). 

4.7  Input  Parameters 

4.7.1  Soil  Properties 

The  soil  used  for  modelling  the  flat  plate-penetration  using  the  MCC  model  is  50-50 
blend  of  kaolinite  and  Edgar  sand.  Figure  4.2  presents  the  grain  size  distribution  of 
both  soils.  It  is  apparent  from  the  distribution  that  the  soil  borders  between  clay 
and  silt.  The  results  of  Atterberg  Limit  tests  (performed  separately  on  Kaolinite  and 
sand  and  also  on  their  mixture  yielded  the  following  results: 

A  series  of  controlled  gradient  consolidation  (CGC)  tests  and  isotropically  consoli¬ 
dated  triaxial  (CIU)  tests  were  performed  to  provide  the  input  parameters  for  MCC 
model.  In  the  CGC  tests,  a  small  pore  pressure  difference  was  maintained  between 
the  two  end  of  the  specimen  by  varying  the  axial  load.  The  soil  sample  was  saturated 
under  a  413kPa  back  pressure.  An  excessive  20.68kPa  pore  presure  difference  was 
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Grain  Diameter,  mm 


Figure  4.2:  Grain  size  distribution  curves 
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maintained  at  the  lower  face  of  the  specimen. 

The  CIU  triaxial  tests  were  perfomed  on  soil  specimens  normally  consolidated  under 
confining  stresses  of  218  kPa,  360  kPa  and  505  kPa.  The  specimens  were  saturated 
under  a  back  pressure  of  414  kPa.  Details  of  the  test  procedures  and  results  of  the 
consolidation  and  triaxial  tests  !x  !xare  presented  in  the  Chapter  on  "Reference  Soil 
Properties”. 

4.7.2  MCC  Model  Parameters 

The  MCC  model  requires  two  consolidation  parameters,  namely  A  and  k  that  are 
related  to  the  conventional  parameters  Cc  and  C,  as: 

k  =  0.434C,  (4.50) 

A  =  0.434Cc  (4.51) 

The  model  also  uses  a  strength  parameter,  M,  the  shear  elastic  modulus,  G,  and  the 
initial  void  ratio  eo- 

The  strength  parameter,  M,  was  evaluated  from  the  stress  path  obtained  from  an 
undrained  triaxial  test,  Figure  4.3.  The  slope  of  the  tangent  to  this  curve,  starting 
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Input  parameters  for  the  Modified  Cam- Clay  model 


PARAMETER 

USED  FOR 

HOW  TO  OBTAIN 

M  or  c 

(M  =  0.93) 

(c  =  0.759) 

Shape  of  Yield  Surface 

M  equals  the  common  tan¬ 
gent 

to  the  stress  paths  of 

undrained 

triaxial  tests  with  different 

preconsolidation  stresses 
c=^2/3M 

Cc  or  A 

(A  =  0.13) 

(Ce  =  .299) 

Changes  of  the  size  of 

the  Yield  Surface 

From  the  virgin  isotropic 

consolidation  curve 

A  =  (0.434).Cc  —  de/d(lrio) 

C „  or  k 

(k  -=  0.026) 

(C,  =  0.0599) 

Changes  of  the  size  of 

the  Yield  Surface 

From  the  isotropic  rebound 

curve: 

k  =  (0.434).CS  =  de/d(lria) 

G 

(=  130  Kpa) 

Elastic  Shear  Modulus 

From  the  small  strain  region 

of  shear  stress-strain  curve 

of  undrained  triaxial  test 

c0(=  0.6889) 

Initial  void  ratio 

Laboratory  Measurement 

Table  4.3:  Input  parameters  for  the  MCC  model 
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from  the  origin  is  equal  to  M. 


The  consolidation  parameters,  A  and  k  were  obtained  from  a  series  of  controlled  gra¬ 
dient  consolidation  tests.  Results  of  these  tests  will  be  presented  in  Chapter  7.  The 
slope  of  the  virgin  compression  curve  equals  A  and  the  average  slope  of  the  unload- 
reload  curve  equals  k. 

The  Elastic  shear  modulus  was  obtained  from  the  shear  stress-strain  curve  Figure  4.4. 
G  equals  the  slope  of  this  curve  in  the  small  strain  region.  The  parameters  used  in 
the  MCC  model  and  their  values,  as  obtained  from  laboratory  tests  are  presented  in 
table  4.3.  The  calibrated  model  was  used  to  generate  the  shear  stress-strain  curve  of  a 
undrained  triaxial(CIU)  test  on  the  soil  (Figure  4.4).  The  MCC  Model  overestimates 
the  peak  shear  stress  (<r(  —  03)  by  16%.  This  discrepancy  is  not  entirely  unexpected 
since  the  test  soil  consists  of  50%  sand  whereas  the  MCC  Model  is  more  suitable  for 
soils  with  higher  clay  content.  Sand  was  used  for  speeding  up  the  drainage  during 
consolidation  of  a  soil  sample  in  the  calibration  chamber.  Thus,  using  the  model 
parameters  described  above,  the  effective  stress  increment  (or  its  components  F 
and  Sij)  caused  by  a  given  strain  increment  e.j  can  be  evaluated. 
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shear  stress,  (a'l-a'a)/!^ 


mean  effective  stress,  (cr^-f <j'2+cr'3)/3g' 


Figure  4.3:  Stress  paths  for  MCC  mode]  and  triaxial  tests. 
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shear  stress,  (oY-oY)/^' 


Figure  4.4:  Shear  stress  versus  shear  stiain  for  MCC  mode]  and  triaxial  test. 
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Chapter  5 

Deviatoric  Stress  Field 


The  strain  path  along  a  streamline  forms  the  input  to  the  MCC  model.  Computations 
resulted  in  the  corresponding  change  of  deviatoric  stress  ( stJ )  and  mean  effective 
stress  ( <7 ).  In  demonstrating  the  stress  field,  the  same  reference  radius  R  described  in 
Chapter  2  is  used  to  denote  the  relative  position  of  the  soil  element.  R  is  defined  as 
half  the  minimum  dimension  of  the  penetrometer.  The  following  R  values  were  used 
for  different  penetrometers: 

Rdmt  =  7mm 
Rcone  —  17.8mm 

R Stepped  Blade  *  1.15mm 

R-Tapered  Blade  3.3m.m 

The  same  coordinate  system  chown  in  Figure  3.2  is  used  to  ptresent  the  results. 

5.0.1  5t  Stress  Space 

In  three  dimensional  problems,  a  stress  state  is  described  by  six  independent  compo¬ 
nents.  Since  the  flat-plate  penetration  deformation  is  posed  as  an  undrained  shearing 
problem,  the  sum  of  normal  stress  components  vanishes  (cr„  =  0).  The  soil  defor¬ 
mation  mechanism  can  not  be  simulated  by  any  one  simple  soil  test.  It  is,  however, 
possible  to  view  the  soil  deformation  as  a  combination  of  shearing  modes  in  Triaxial 
Compression  test  (TC),  Direct  Simple  Shear  (DSS)  and  Pressuremeter  Test  (PMT) 
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(Prevost,  1978).  To  relate  the  soil  deformation  mechanism  to  the  above  mentioned  soil 
tests  an  Si  subspace  similar  to  the  one  suggested  by  Kavvadas  (1982)  is  adopted.  The 
components  of  5-dimensional  5,  stress  space  are  related  to  components  of  deviatoric 
stress  tensor  as  follows: 

51  =  (2 sw  -  sxx  -  szz)/V6 

52  =  (Szz~Sxx)/V2 

53  =  Syz/y/2 

54  =  sxy/y/ 2 

55  =  Wv/2  (5.1) 

Si  denotes  the,  shearing  mode  in  triaxial  compression,  S 2  signifies  the  cylindrical 
cavity  expansion  mode  (the  cavity  being  infinitely  long  in  the  x  direction)  and  S3 

implies  the  shearing  mode  in  direct  simple  shear.  S4  and  S$  are  similar  to  the  S3  but 

in  the  xy  and  zx  directions  respectively. 

5.0.2  Stress  Distribution 

The  stress  distributions  due  to  different  steady  penetrating  penetrometes  were  ana¬ 
lyzed  and  will  be  presented  in  the  following  text  in  terms  of: 

1.  stress  profiles:  plots  showing  stress  components  versus  y/R; 

2.  stress  paths:  plots  of  5,  versus  S},i  ^  j\ 

3.  stress  contours:  plots  showing  distribution  of  different  stress  components  in  a 
plane. 


5.1  Axisymmetric  versus  Three-dimensional  Pen¬ 
etration 

To  evaluate  the  effect  of  geometry  on  strain  rates,  shear  strain  rate  (%<.*)  profiles 
corresponding  to  z/R  =  1  were  compared  for  the  cone  and  DMT  (Figures  5.1  and 
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5.2).  In  the  case  of  the  DMT,  profiles  were  plotted  for  x/R  values  of  0,  4,  5  and  6. 
The  x/R  value  of  6  corresponds  to  approximately  the  corner  of  the  DMT. 

For  both  penetrometers  strain  rates  are  very  high  (maximum  strain  rate  ss  20%  /sec.) 
in  the  tip  region  (OR  —  5.67i?).  However  there  is  no  further  similarity  between  the 
cone  and  the  DMT.  Shear  strain  rate  shows  a  positive  peak  at  the  base  of  cone  tip 
(y/R  =  5.67),  whereas  7^  profiles  for  the  DMT  are  characterized  by  a  positive  peak 
at  the  tip  (y/R  =  0)  and  a  negative  peak  at  the  base.  Such  a  reduction  in  strain 
rates  to  the  extent  of  forming  a  negative  peak  is  unique  to  flat  plate  geometry.  It 
is  caused  presumably  due  to  the  presence  of  a  limited  blade  thickness  which  allows 
strain  dissipation  or  relaxation  at  the  corners  of  the  blade.  The  influence  of  DMT 
corners  resulting  in  strain  rate  relaxation  around  the  edges  will  hereafter  be  referred 
to  as  the  end  effect.  Shear  strain  rates  for  streamlines  further  from  the  center  of 
the  blade  show  less  reduction  in  the  tip  region.  Such  a  compensation  in  strain  rate 
occurs  because  soil  deformation  in  x  direction,  on  account  of  the  end  effect,  is  more 
significant  towards  the  sides  of  the  blade  (Figure  5.1). 

The  ensuing  text  deals  with  differences  in  the  stress  fields  due  to  axisymmetric  and 
flat  plate  penetration.  The  axisymmetric  character  of  the  cone  makes  it  sufficient 
to  present  stress  in  any  one  vertical  plane  passing  through  its  vertical  axis.  However 
there  is  no  such  representative  plane  for  the  DMT.  Therefore  the  stresses  due  to  DMT 
penetration  were  examined  in  three  xz  planes  and  one  yz  plane. 

1.  y/R  =  0:  horizontal  plane  at  the  tip; 

2.  y/R  =  5.67:  horizontal  plane  at  the  base  of  the  tip; 

3.  y/R  =13:  horizontal  plane  at  the  center  of  the  DMT  membrane; 

4.  x/R  =  0:  vertical  plane  through  the  center  of  the  DMT. 

In  the  MCC  approach  the  DMT  and  cone  penetration  mechanism  are  both  associated 
with  symmetric  deviatoric  stress  tensors  with  four  significant  components,  namely; 
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sXx )  -Syvt  szz  and  syz.  Figure  5.3  shows  the  location  of  the  data  labels  (relative  to  the 
penetrometer)  used  in  marking  stress  path  plots  for  the  cone  and  DMT. 

In  addition  to  the  deviatoric  stress  components,  shear  octahedral  stress  (r^)  and 
mean  effective  stress  a  were  also  used  to  study  the  stress  field  due  to  both  penetrom¬ 
eters.  Plots  are  presented  to  compare  stress  fields  around  a  cone  and  a  DMT.  The 
apex  angle  for  both  penetrometers  being  20°,  the  tapered  region  or  the  tip  extends 
from  0i2  —  5.67 R.  Stress  components  of  both  the  S,  space  and  the  deviatoric  stress 
tensor  were  plotted  in  the  tip  region  of  the  cone  and  the  DMT  (Figures  5.4,  5.5,  5.6 
and  5.7).  The  following  observations  were  made: 

•  sxx,  the  least  significant  of  the  stress  components  for  DMT,  is  the  predominant 
tensile  stress  component  for  the  cone  i.e.  it  reaches  the  maximum  magnitude 
in  tension.  The  smaller  dimension  of  the  DMT  allows  relaxation  of  sxx  at  the 
corners. 

•  Syy  for  both  penetrometers  is  in  compression  just  before  the  tip.  of  DMT 
undergoes  a  complete  unload-reload  cycle  in  the  tip  region.  sw  of  cone  shows 
similar  cyclic  characteristics  but  the  reloading  is  only  partial.  The  vertical  stress 
components  iw  is  pronounced  for  both  penetrometers  but  picks  up  magnitude 
earlier  for  DMT. 

•  szz  shows  an  loading-unloading  cycle  for  both  penetrometers  but  unloading  is 
partial  for  the  cone. 

•  syz  of  cone  and  syz  of  DMT  are  comparable  in  both  magnitude  and  direction. 

The  above  stress  characteristics  can  also  be  observed  and  explained  using  the  5, 
components  in  terms  of  stress  states  typical  to  conventional  laboratory  tests.  Si, 
representing  the  triaxial  state  of  stress  follows  the  same  trend  as  sw  for  both  pen¬ 
etrometers.  The  pressuremeter  mode  of  stress  represented  by  S2  is  more  prominent 
for  the  cone  than  the  DMT  (Figure  5.6).  The  S2  versus  S3  plot  shows  that  the  part 
of  the  curve  representing  the  tip  region  is  almost  flat  for  the  DMT.  S3,  which  is  es¬ 
sentially  similar  to  syz  does  not  show  marked  differences  between  the  cone  and  the 
DMT. 
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5.2  Effects  of  Flat  Plate  Geometry 


szz  and  Tod  contours  for  DMT  penetration  (Figures  5.10  and  5.11)  show  that  towards 
the  end  of  the  blade  (in  x  direction),  a  sharp  stress  relaxation  occurs.  This  phenomena 
is  obviously  not  accounted  for  by  plane  strain  approximation  (i.e.  w/t  =  oo).  Such 
a  redistribution  of  stress  around  corners  is  expected  for  all  fiat  plates  and  is  due 
to  the  end  effect.  The  extent  and  influence  of  end  effect  and  its  dependence  on 
shape  geometry  of  the  penetrometer  was  examined.  Analyses  were  carried  out  for  the 
following  variations: 

1.  Different  w/t  ratios,  where  w  is  the  width  and  t  is  the  thickness  of  the  flat 
plate.  In  all  cases  t  was  equal  to  14mm  and  the  angle  of  taper  was  20°  which 
are  the  same  as  DMT.  These  flat  plates  will  be  referred  to  as  the  simple  flat 
plates  herein  after; 

2.  Other  flat-penetrometers,  namely  a)  Tapered  Blade  and  b)  Stepped  Blade. 

To  illustrate  the  effect  of  lateral  distance  ( zj R )  on  a  stress  element,  stress  paths  are 
plotted  for  elements  at  1J?  and  5jR.  A  comparison  of  stress  paths  for  these  two  el¬ 
ements  shows  that  points  further  away  experience  a  smaller  range  of  stress  states. 
Figures  5.12  shows  that  differences  between  a  three  dimensional  simple  flat  plate  and 
its  two-dimensional  (plane  strain  or  w/t  =  oo)  approximation  can  be  reduced  by 
increasing  w/t  ratio,  w/t  =  20  was  found  to  have  stress  paths  very  similar  to  two- 
dimensional  stress  paths.  Increasing  the  w/t  ratio  makes  the  center  relatively  closer 
to  the  side  of  the  blade  thereby  increasing  the  influence  of  end  effect  on  the  stress 
path.  Stress  paths  for  lower  w/t  ratios  compare  closer  with  those  caused  by  cone 
penetration  as  opposed  to  the  plane  strain  approximation. 

To  further  demonstrate  the  three-dimensional  nature  of  the  flat  plate  penetration, 
stress  paths  and  stress  profiles  of  three  different  penetrometers  (DMT,  Tapered  Blade 
and  Stepped  Blade)  were  plotted  (Figures  5.14  -  5.21  ).  The  stress  paths  for  tapered 
blade  are  essentially  similar  to  those  for  the  DMT  because  of  the  inherent  similarity 
in  shape.  Stress  paths  and  stress  profiles  for  stepped  blade  show  a  cyclic  character 
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because  of  the  loading-unloading  of  stresses  caused  by  steps.  The  penetration  mech¬ 
anism  of  the  stepped  blade  can  be  thought  of  as  a  continuous  penetration  of  a  series 
of  flat  plates  (of  the  DMT  geometry)  with  decreasing  w/t  ratios.  The  Tapered  Blade, 
according  to  design  by  Handy  and  Lutenegger  (1985)  has  a  range  of  w/t  ratios  from 
13.33  at  the  tip  to  7.33  at  the  end.  The  Stepped  Blade  (Handy  et  al.,  1992),  has  four 
steps  with  w/t  ratios  of  10.8,  7.14,  5.4  and  4.31  respectively.  Both  Stepped  Blade 
and  Tapered  Blade  have  high  initial  w/t  ratios  which  explains  the  similarity  between 
the  initial  parts  of  2D  and  three  dimensional  stress  paths. 


76 


shear  strain  rate,  (7oct) 


y/R 


Figure  5.1:  Shear  strain  rate  (7oc£/sec)  profiles  for  cone  and  DMT. 
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Figure  5.4:  Profiles  of  normalized  stress  components  for  cone  penetration. 
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Figure  5.6:  Normalized  stress  paths  for  20°  cone. 
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Figure  5.8:  szx  contours  for  cone  and  DMT. 
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Figure  5.9:  roc(  contours  for  cone  and  DMT. 
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Figure  5.10:  szz  contours  for  DMT  in  different  horizontal  planes. 
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Figure  5.12:  Stress  paths  showing  effect  of  t eft  ratio. 
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Figure  5.13:  Stress  paths  for  different  x/R  values. 
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Figure  5.16:  Stress  profiles 


a  stepped  blade  penetration. 
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Figure  5.17:  Stress  paths  for  two-dimensional  DMT 
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Figure  5.18:  Stress  paths  for  two-dimensional  tapered  blade 


Figure  5.19:  Stress  paths  for  three-dimensional  tapered  blade. 
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Figure  5.20:  Stress  paths  for  two-dimensional  stepped  blade. 


96 


Figure  5.21:  Stress  Paths  for  three-dimensional  tepped  blade. 
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Chapter  6 

Total  Stress  and  Pore  Pressure 


Since  the  soil  was  treated  as  incompressible,  normal  octahedral  stress  aoct  can  not 
be  estimated  from  MCC  model  and  therefore  must  be  determined  from  equilibrium 
considerations.  Asuming  the  clay  as  fully  saturated,  the  mean  normal  stress  related 
pore  pressure  varations  can  be  estimated  from  aoct.  Along  with  the  deviatoric  stresses, 
estimates  of  o0Ct  and  pore  pressure  comprise  of  a  complete  solution  to  the  penetration 
problem. 


6.1  Equilibrium  Considerations 

For  a  general  three-dimensional  penetration  problem  the  equilibrium  in  x,  y  and  z 
directions  can  be  written  as: 


d&oct 

()Sxx 

dSxy 

d$X2 

dx 

dx 

dy 

dz 

1 

Tj 

o 

b 

<0 

OSyx 

dsyy 

dSyz 

dy 

dx 

dy 

dz 

9(7  oct 

ds:x 

dS-y 

dszz 

dz 

dx 

dy 

dz 

For  the  axisymmetric  penetration  problem  the  above  set  equations  can  be  reduced  in 
the  cylindrical  coordinate  system  as: 

dcroct  _  cbs  rr  9srz  Srr  $00 

dr  dr  dz  r 
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3<Joct 

~dT 


3srz  3szz  $r2 
dr  8  z  r 


(0.2) 


where  the  r  in  cylindrical  coordinates  corresponds  to  2  direction  in  cartesian  system. 
9  in  cylindrical  coordinates  corresponds  to  x  direction  in  cartesian  system  and  2  in 
cylindrical  coordinates  corresponds  to  y  direction  in  cartesian  system. 


Levadoux  and  Baligh  (19S0)  showed  that  integration  of  any  one  of  the  set  of  equations 
6.2  yields  the  value  of  aoct.  The  value  of  aoct  thus  determined  is  path  dependent  i.e. 
it  is  different  for  integration  in  r  and  2  direction,  unless  the  following  condition  is 


satisfied. 


3  .  3aoct  8  daoct 

8r[  8z  ’  ~  Dz  dr  ’ 


(6.3) 


The  above  condition  can  be  satisfied  only  if  the  strains  are  compatible  with  the  devi- 


atoric  stress  model.  Furthermore,  strains  should  be  determined  as  an  exact  solution 


to  the  penetration  problem  rather  than  approximate  solutions  corresponding  to  an 
ideal  fluid. 


The  simplistic  extension  of  the  concept  used  by  Baligh  (1985),  to  the  general  three- 
dimensional  case  of  DMT  penetration  proved  to  be  less  desirable.  The  error  in  in¬ 
tegration  is  multiplied,  presumably  because  of  three  directions  in  which  integration 
can  be  performed  to  calculate  aoct.  To  correctly  analyze  this  problem  it  is  neces¬ 
sary  to  solve  a  poisson  equation  in  three-dimension.  Such  computation  is  beyond  the 
scope  of  this  research.  The  following  evaluation  of  pore  pressure  and  total  stresses 
are  therefore  restricted  to  plane  strain  cases.  Since  the  equilibrium  considerations 
have  been  successfully  employed  to  determine  aoct  for  axisymmetric  cone  penetration 
(Teh,  1987;  Baligh,  19S5),  it  is  therefore  expected  that  two-dimensional  plane  strain 
approximations  of  the  flat  penetrometers  will  generate  reasonable  estimates  of  croct 
and  pore  pressure.  In  plane  strain  simplification,  the  width  of  the  penetrometer  is 
assumed  to  be  infinite  in  x  direction  (Figure  3.2).  The  equilibrium  equations  for  plane 
strain  simplification  can  be  written  as: 

8(XoCt  8  s  yy  3SyZ 

3y  3y  dz 
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d(T0ct 

dz 


(6.4) 


Oszy  Q  &  zz 

~di~~dT 

where  integration  can  be  performed  in  the  y  and  2  direction  (Figure  3.2).  A  centra] 
difference  iterative  integration  scheme  was  used  to  evaluate  eroc(.  The  scheme  utilized 
isoparametric  elements  (Levadoux  and  Baligh,  1980)  for  evaluating  derivatives.  The 
fluid  flow  idealization  of  soil  behavior  during  deep  steady  penetration  is  believed 
generally  to  overestimate  the  shear  strains  ey2  (or  er2  in  cylindrical  coordinates).  This 
is  because  of  the  lack  of  shearing  resistance  in  a  perfect  fluid.  Furthermore,  in  the 
deviatoric  stress  model  the  predicted  shear  stress  syz  (or  sr2  in  cylindrical  coordinates) 
is  controlled  by  shear  strain  £yz  (or  eTZ  in  cylindrical  coordinates)  which  can  therefore 
be  unrealistic.  It  was  shown  by  Baligh  (1985)  and  Teh  (1987)  that  for  integration 
in  2  direction,  the  effect  of  syz  is  minimized  and  hence  more  reliable  estimates  of 
cr0Ct  are  obtained.  Therefore,  prediction  of  pore  pressure  and  total  stresses  presented 
subsequently  are  determined  by  means  of  integration  along  the  2  direction. 


6.2  Results 

Pore  pressure  and  total  stress  computations  were  performed  for  the  various  penetrom¬ 
eters  discussed  so  far.  In  addition,  a  flat-plate  geometry  similar  to  a  two-dimensional 
DMT  but  with  a  60°  tip  was  also  analyzed  in  order  to  study  the  effect  of  varying  the 
tip  angle. 

Figures  6.1  and  6.2  show  the  pore  pressure  and  total  stress  profiles  for  a  20°  DMT 
and  a  cone  respectively.  It  is  evident  from  these  figures  that  for  the  same  apex  angle 
the  total  stress  state  just  ahead  of  the  tip  is  different  for  axisymmetric  (cone)  and 
plane  strain  (two-dimensional  DMT)  penetration.  Pore  pressure  and  total  stress  is 
generally  higher  for  cone  penetration.  Furthermore,  the  rate  of  change  of  total  stress 
and  pore  pressure  is  higher  for  axisymmetric  penetration.  Since  the  end  effect  in 
the  case  of  three-dimensional  DMT  causes  stress  relaxation  around  corners  of  the 
flat-plate,  peak  pore  pressure  and  total  stresses  can  be  expected  to  be  lower  than  the 
plane  strain  approximation.  The  pore  pressure  due  to  axisymmetric  and  flat-plate 
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penetration  being  different  the  initial  conditions  for  the  subsequent  consolidation  are 
also  expected  to  be  dissimilar. 

The  effect  of  shape  and  geometry  was  studied  by  analyzing  flat-plate  geometries  with: 

1.  a  60°  tip; 

2.  a  3°  taper  (Tapered  Blade);  and 

3.  steps  (Stepped  Blade) 

Pore  pressure  and  total  stress  computations  for  the  above  mentioned  flat  plate  ge¬ 
ometries  were  carried  out.  The  results  have  been  presented  at  the  end  of  this  chapter 
in  form  of  profiles  an  contours  of  pore  pressure  and  total  stresses.  Characteristics  of 
pore  pressure  and  total  stress  distribution  can  be  summarized  as  follows: 

•  Pore  pressure  and  total  stress  profiles  (Figures  6.1  -  6.5)  show  a  positive  peak 
in  the  tip  region  for  all  penetrometer  geometries. 

•  The  peak  magnitude  for  pore  pressure  and  total  stress  is  comparable  for  the 
DMT  and  the  cone. 

•  Increasing  the  tip  angle  to  60°  for  the  DMT  makes  the  total  stresses  more 
concentrated  in  the  tip  region,  however  there  is  no  significant  increase  in  the 
peak  magnitude.  Furthermore,  the  total  stresses  and  pore  preessure  fall  quicker 
(after  the  tip)  and  to  a  lower  value  when  the  tip  angle  was  increased  to  60°. 

•  As  observed  in  the  case  of  deviatoric  stress  distribution,  pore  pressure  and 
total  stresses  for  tapered  blade  and  stepped  blade  do  not  exibit  an  exponential 
increase  with  the  thickness  of  the  blade.  The  total  stress  and  pore  pressure 
profile  for  the  tapered  blade  shows  a  sharp  peak  of  magnitudes  significantly 
higher  than  those  observed  for  the  DMT  and  cone.  Profiles  for  the  stepped 
blade  show  peaks  corresponding  to  the  loading  unloading  cycles  caused  due  to 
each  step  of  increasing  w/t  ratio.  Pore  pressure  and  the  different  components 
of  total  stress  show  almost  identical  trends  for  both,  stepped  blade  and  tapered 
blade. 
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Figure  6.1:  Normalized  total  stress  and  pore  pressure  profiles  for  20°  DMT  . 
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Figure  6.2:  Normalized  total  stress  and  pore  pressure  profiles  for  20°  cone  . 
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Figure  G . 4 :  Normalized  total  stress  and  pore  pressure  profiles  for  tapered  blade. 
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Figure  6.5:  Normalized  total  stress  and  pore  pressure  profiles  for  stepped  blade 
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Figure  6.7:  A u/P'0  contours  for  20°  DMT. 
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Figure  6.9:  crzzj  P'Q  contours  for  20°  DMT. 
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Figure  6.11:  A u/P'0  contours  for  20° 
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cone. 
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Figure  6.13:  azzj P'0  contours  for  20°  DMT. 
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Figure  6.17:  azz/ P'Q  contours  for  60°  DMT. 
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Figure  6.19:  Au/P°  contours  for  tapered  blade. 
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Figure  6.21:  azJP "0  contours  for  tapered  blade. 
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Figure  6.23:  A ujP'0  contours  for  stepped  blade. 


115 


2  /R 

Figure  6.24:  <ryy/Pi '0  contours  for  stepped  blade. 
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Figure  6.25:  <yzzjP'0  contours  for  stepped  blade. 
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Chapter  7 


Reference  Soil  Engineering 
Properties 


A  mixture  of  509c  kaolinite  and  50%  Edgar  sand  was  used  to  create  t lie  clay  specimens 
in  the  calibration  chamber  tests.  To  establish  the  geotechnical  engineering  properties 
of  the  clay,  a  series  of  laboratory  experiments  were  performed.  In  addition  to  the  de¬ 
termination  of  basic  engineering  properties,  the  experiments  included  consolidation 
and  triaxial  tests.  This  chapter  describes  the  details  of  those  experiments  and  their 
results. 


7.1  Soil  Classification 

Figure  7.1  shows  the  grain  size  distributions  of  the  kaolinite  and  Edgar  sand.  At- 
terberg  limit  tests  (ASTM  D 4 3 1 S )  were  conducted  on  the  1 00 9c  kaolinite  and  the 
kaolinite/Edgar  sand  mixture  with  the  results  listed  in  Table  7.1.  In  this  table  the 
specific  gravity  (ASTM  DS50)  of  kaolinite  and  Edgar  sand  is  also  shown.  The  100% 
kaolinite  will  be  referred  to  as  K100  and  the  50/50  mixture  will  be  referred  to  as  K50 
hereinafter. 

Based  on  1JSCS  the  kaolinite  can  be  classified  as  MII-C1I  which  is  typical  for  kaoli- 
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Figure  7.1:  Grain  size  distribution  of  kaolinite  and  Edgar  sand. 


Soil 

Type 

Liquid 

Limit  (%) 

Plastic 

Limit  (%) 

Plastic 

Index  (Vc) 

Specific 

Gravity 

Kaolinite 

51.5  to 

25.2  to 

23.4  to 

2.62  to 

56.0 

30.1 

26.3 

2.69 

Kaolinite  and 

29.2  to 

16.3  to 

12.5  to 

N/A 

Sand  (50/50  Mix) 

30.9 

16. S 

14.1 

Sand 

N/A 

N/A 

N/A 

2.63  to 

2.71 

Table  7.1:  Engineering  properties  of  K50  soil. 


nite  (Mitchell,  197G).  The  Edgar  sand  was  classified  as  SP  or  poorly  graded  fine  sand. 


7.2  Consolidation  Characteristics 

To  provide  the  consolidation  characteristics  of  the  clay  specimens  used  in  the  re¬ 
search,  a  series  of  controlled  gradient  consolidation  (CGC)  tests  were  conducted.  The 
advantages  of  CGC  method  are  well  recognized  (Lowe  et  ah,  1969).  This  method 
is  especially  attractive  since  it  can  be  servo-controlled  by  a  computer  and  therefore 
eliminates  the  need  to  take  readings  and  to  adjuste  the  consolidation  pressure  man¬ 
ually. 


7.2.1  The  Controlled  Gradient  Consolidation  Test  Proce¬ 
dures 

The  clay  specimens  were  prepared  in  the  laboratory  from  slurry  following  the  tech¬ 
niques  developed  by  Huang  et  al.  (19SS).  The  slurry  was  mixed  with  an  initial 
water  content  of  twice  the  liquid  limit.  Deionized  and  deaired  water  was  used  in  the 
preparation  of  the  slurry.  Upon  mixture,  the  slurry  was  consolidated  under  a  vertical 
pressure  of  138  kPa.  At  the  end  of  slurry  consolidation,  a  clay  specimen  with  suffi¬ 
cient  consistancy  to  support  its  own  weight  was  created. 

The  controlled  gradient  consolidation  testing  device  fabricated  as  part  of  the  re¬ 
search  is  shown  in  Figure  7.2.  The  specimen  is  completely  contained  to  enable  a  back 
pressure  to  be  applied  to  the  specimen.  In  all  the  CGC  tests,  a  back  pressure  of  414 
kPa  was  applied.  A  Bdlofram/Ymear  ball  bearing  piston  configuration  was  used  to 
minimize  shaft  friction.  The  soil  specimen  was  obtained  by  pushing  a  tapered  ring 
(angle  of  taper  =  3°)  into  the  clay.  The  ring  has  an  inside  diameter  of  61  mm  and  a 
height  of  17.8  mm. 
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Figure  7.2:  Cross  sectional  view  of  the  controlled  gradient  consolidation  device. 
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The  consolidation  pressure  was  monitored  by  pressure  tranducers  located  at  the 
top  and  bottom  of  the  specimen.  This  arrangement  enables  any  significant  friction  on 
the  side  of  the  specimen  to  be  detected.  Drainage  was  allowed  only  through  the  top  of 
the  specimen  against  the  back  pressure.  A  differential  pressure  transducer  was  used 
to  measure  the  pore  pressure  difference  between  the  top  and  bottom  of  the  specimen. 
The  vertical  consolidation  pressure  was  servo-controlled  by  a  computer  to  maintain  a 
positive  excess  pore  pressure  of  20.7  kPa  at  the  bottom  of  the  specimen.  A  Fairchild 
T5221  I/P  and  E/P  transducer  regulator  was  used  to  control  the  consolidation  pres¬ 
sure. 


7.2.2  Consolidation  Test  Results 

Plots  of  the  test  data  are  shown  in  Figures  7.3  to  7.10.  In  tests  K50-2  and  K100-2 
the  consolidation  pressure  was  increased  to  276  kPa  and  then  unloaded  prior  to  the 
controlled  gradient  consolidation  tests.  The  preconsolidation  pressures  depicted  in 
either  the  lime-displacement  (Figures  7.3  and  7.7)  or  the  e-log  p  curves  (Figures  7.6 
and  7.10)  are  essentially  the  same  as  the  known  values  listed  in  Table  7.2.  This  is  an 
indication  of  the  good  quality  of  the  tests.  Table  7.2  reveals  key  parameters  found  in 
the  CGC  tests.  The  Cv  values  are  obtained  by  taking  derivatives  of  laboratory  data. 
The  somewhat  erratic  results  are  therefore  expected  and  numbers  given  in  Table  7.2 
represent  a  range  of  Cv  values  for  the  indiviual  tests.  The  K50  has  a  Cv  typically  at 
2  to  5  times  that  of  K100.  This  is  clearly  due  to  the  presence  of  sand  in  the  K50 
specimens. 
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Preconsolid. 

10  3cm2/sec 

Test  No. 

stress,  kPa 

Virgin  Loading 

Unload 

Reload 

I\  100-1 

13S 

1.5  to  3.5 

2.0  to  2S.0 

5.0  to  15.0 

K 100-2 

276 

1.5  to  3.5 

2.0  to  28.0 

5.0  to  15.0 

K50-1 

13S 

2.0  to  7.0 

10.0  to  50.0 

5.0  to  40.0 

K50-2 

276 

2.0  to  7.0 

10.0  to  50.0 

5.0  to  40.0 

Table  7.2:  Summary  of  the  controlled  gradient  consolidation  tests. 
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Figure  7.5:  C„  versus  effective  consolidation  stress  for  K 100-2. 
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Figure  7.9:  Cv  versus  effective  consolidation  stress  for  K50-2. 


7.3  Undrained  stress-strain  and  Strength  Rela¬ 
tionships 

A  series  of  isotropically  and  anisotropically  consolidated  undrained  triaxial  tests  (CIU 
and  CAU)  were  conducted  to  establish  the  stress-strain  and  strength  relationships  of 
the  K50  clay.  Triaxial  tests  were  limited  to  K50  specimens  since  this  was  the  only 
soil  used  in  the  calibration  chamber  tests. 

The  same  slurry  preparation  and  consolidation  process  as  in  the  consolidation 
tests  were  used  to  create  triaxial  test  specimens.  The  slurry  consolidometer  was 
made  of  a  plexiglass  tube  with  an  inside  diameter  of  62  mm  and  a  total  height  of  380 
mm  (see  Figure  7.11).  Consolidation  pressure  was  applied  through  the  top  piston. 
Drainage  was  allowed  at  both  the  top  and  bottom  of  the  specimen.  The  height  of  the 
specimens  were  approximately  150  mm  and  the  water  content  was  at  27-28%  at  the 
end  of  slurry  consolidation.  Upon  extrusion,  the  specimen  was  inserted  directly  into 
a  rubber  membrane.  Trimming  was  necessary  only  at  the  ends  of  the  specimen  to  a 
final  required  height  for  the  triaxial  test  of  140  mm.  All  specimens  were  consolidated 
and  sheared  under  a  back  pressure  of  550  kPa. 


7.3.1  CIU  Triaxial  Test  Results 

The  CIU  tests  included  axial  compression  and  extension  tests  on  specimens  normally 
consolidated  ( a'v0  =  218,  360,  505  kPa)  and  over  con  soli  da  ted  (OCR  =  5  and  10).  To 
evaluate  the  strain  rate  effects,  a  series  of  tests  were  performed  on  normally  consol¬ 
idated  specimens  with  strain  rates  of  0.013,  0.127,  and  0.508  cm/min.  Tests  with 
unload  reload  cycles  were  performed  on  NC  and  OC  specimens.  These  cycles  were 
applied  after  the  peak  stress  difference  was  reached.  The  SPM  studies  indicated  that 
soil  elements  around  a  penetrometer  typically  experience  a  process  of  unloading  (cone 
and  flat-plate)  and  reloading  (flat-plate)  only  after  a  large  strain  or  soil  failure  has  oc¬ 
curred.  Previous  laboratory  experimental  studies  have  been  limited  to  unload-reload 
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cycles  prior  to  soil  failure.  Therefore,  these  cyclic  tests  have  important  implications 
for  the  interpretation  of  penetration  tests. 


A  number  of  key  parameters  were  determined  for  each  triaxial  test.  They  in¬ 
cluded  Skempton’s  pore  pressure  parameter  at  failure  (A/)  (Skempton,  1954),  and 
the  undrained  Young's  modulus  at  50  %  of  the  maximum  deviator  stress  (£50).  Other 
key  parameters  included  su,  the  undrained  shear  strength,  <f>',  the  friction  angle  and 
M,  the  slope  of  the  critical  state  line.  These  parameters  are  defined  as  follows. 


Af  = 


An 


Eso  = 


Su  = 


A(<7i  -  <T3) 

~  <73)max/2 

Sf/2 

(^1  ^3)mai. 


<t>'  =  sin  1(tan(rJ;)) 


M  =  3(irSy)  +  3  _  (CTi  -  fjlgifal 
2(r^>  +  4  W  +  2»4)/3 


(7.1) 

(7.2) 

(7.3) 

(7.4) 

(7.5) 


Figure  7.12  shows  the  effective  stress  paths  for  the  tests  on  NC  and  OC  speci¬ 
mens.  The  stress-strain  curves  normalized  with  respect  to  the  preconsolidation  stress 
for  all  the  tests  on  NC  specimens  are  shown  in  Figure  7.13.  The  closeness  of  these 
normalized  stress-strain  curves  indicates  that  the  tests  are  reasonably  repeatable  and 
that  K50  has  a  fairly  normalized  behavior.  Tables  7.3  and  7.4  present  a  summary 
of  the  key  parameters  found  in  compression  and  extension  triaxial  tests  respectively. 
The  ( A j)  values  would  indicate  that  the  K50  may  behave  like  a  sensitive  clay.  This  is 
most  likely  due  to  the  substantial  amount  of  sand  in  the  clay  mixture.  However,  the 
stress-strain  curves  show  only  slight  strain  softening  in  the  axial  compression  tests  on 
NC  specimens  (see  Figure  7.13). 


Figures  7.14  and  7.15  show  the  effective  stress  path  of  the  tests  with  load-unload 
cycles  on  NC  and  OC  specimens  respectively.  The  data  show  that  regardless  of  the 
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stress  history,  the  excess  pore  pressure  experiences  an  increase  and  then  decrease  dur¬ 
ing  the  reloading  process  (Figures  7.16  and  7.17).  The  stress  path  during  reloading  is 
similar  to  that  of  an  OC  clay  even  when  the  soil  is  normally  consolidated  as  shown  in 
Figure  7.14.  As  a  result,  the  stress  path  crosses  the  critical  state  line  which  violates 
one  of  the  fundamental  concepts  in  critical  state  soil  mechanics.  It  should  be  empha¬ 
sized  at  this  point  that  the  triaxial  test  is  known  for  its  unreliability  while  performing 
tests  with  unload-reload  cycles.  This  is  mainly  due  to  the  inevitable  delay  during  the 
change  of  loading  direction.  These  tests  may  therefore  not  be  repeatable.  However, 
the  evidence  is  sufficient  to  indicate  that  a  post  failure  unload-reload  process  does 
involve  significant  variation  of  pore  pressure.  Furthermore,  this  variation  may  not 
follow  the  conventional  pattern  as  indicated  by  previous  studies  which  are  limited  to 
unload-reload  cycles  applied  before  failure.  Since  excess  pore  pressure  is  an  impor¬ 
tant  part  of  the  interpretation  and  essentially  all  penetration  tests  involve  post  failure 
unloading  and/or  reloading.  It  is  therefore  imperative  to  establish,  experimentally, 
the  behavior  of  soils  under  these  conditions  before  the  theories  can  be  refined. 

Figure  7.18  shows  the  normalized  stres-strain  curves  of  tests  on  NC  specimens 
with  varying  strain  rates.  The  results  show  an  approximate  increase  of  peak  deviator 
stress  of  17.5%  for  every  ten  fold  increase  of  the  strain  rate. 


7.3.2  CAU  Triaxial  Test  Results 

The  CAU  tests  included  an  axial  compression  and  extension  test  on  specimens  nor¬ 
mally  consolidated  at  218  kPa.  Anisotropic  consolidation  was  performed  while  main¬ 
taining  a  K  value  of  0.52  (<r'h/<r'v).  The  K  value  was  determined  from  the  isotropically 
consolidated  triaxial  tests  by  using  the  4>'  value  and  the  equation,  K  =  1  —  sinfi. 
Drainage  was  allowed  against  the  back  pressure  at  the  top  of  the  specimen  only.  The 
cell  and  axial  pressure  were  servo-controlled  to  maitain  an  excess  pore  pressure  of  34 
kPa  at  the  bottom  of  the  specimen  while  keeping  the  K  value  at  0.52.  Undrained 
shearing  was  performed  at  a  strain  rate  of  0.01%  per  minute.  Key  parameters  deter- 
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mined  from  the  tests  are  summarized  in  Table  7.5.  Plots  of  deviator  stress  and  pore 
pressure  versus  the  axial  strain  are  shown  in  Figures  7.19  and  7.20.  The  strain  paths 
for  both  tests  are  shown  in  7.21. 
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Test 

E9 

Esq,  kPa 

su,  kPa 

Cn 

O 

Co 

c 

sJP'o 

NC(218  kPa) 

1.25 

48,000 

60 

800 

0.28 

NC(360  kPa) 

1.44 

91,000 

91 

1000 

0.25 

NC(505  kPa) 

1.36 

87,900 

123 

700 

0.24 

OC(OCR  =  5) 

0.12 

480,000 

96 

500 

0.89 

OC(OCR  =  10) 

-0.19 

6,600 

64 

100 

1.35 

<f>'  =  29.2°,  M  =  1.194 

Table  7.3:  Summary  of  the  CIU  triaxial  tests  sheared  in  axial  compression. 


Test 

eh 

ESo ,  kPa 

su,  kPa 

Eso/Su 

s»IP'o 

NC(218  kPa) 

0.41 

65,700 

46 

1400 

0.21 

NC(360  kPa) 

0.50 

118,600 

83 

1400 

0.24 

NC(505  kPa) 

0.53 

151,200 

108 

1400 

0.21 

OC(OCR  =  5) 

-0.10 

106,000 

53 

2000 

0.65 

OC(OCR  =  10) 

0.40 

73,000 

73 

1000 

1.04 

(p'  =  28.5°,  M  =  1.189 

Table  7.4:  Summary  of  the  CIU  triaxial  tests  sheared  in  axial  extension. 


A f  E$o,  kPa  su,  kPa  E5o/su 
0.59  110,500  65  1700 

1.11  15,000  20  750 

<f>'  =  27.2°,  M  =  1.180 

Table  7.5:  Summary  of  the  CAU  triaxial  tests  (NC  218  kPa). 
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Figure  7.15:  Effective  stress  path  of  the  OC  load  -  unload  triaxial  test. 
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Figure  7.16:  Change  of  pore  pressure  versus  strain  for  the  NC  load-unload  triaxial 
test. 
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Figure  7.17:  Change  of  pore  pressure  versus  strain  for  OC  load-unload  triaxial  test. 
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Chapter  8 


The  Calibration  Chamber  System, 
Testing  Procedure  and  Program 


A  slurry  consolidation  and  a  calibration  chamber  system  were  fabricated  for  the  re¬ 
search.  The  system  enables  a  clay  specimen  to  be  consolidated  from  slurry  and  then 
transported  to  the  chamber  for  the  second  stage  consolidation  and  the  calibration 
test.  Designs  of  the  testing  devices  and  set  up  followed  the  techniques  developed  by 
Huang  (1986).  The  system  included  a  slurry  consolidometer,  a  computer  servo  con¬ 
trolled  calibration  chamber  system  and  model  fiat-plate  penetrometers.  This  chapter 
provides  the  details  of  these  apparatuses  their  applications  and  the  testing  program. 


8.1  The  Slurry  Consolidometer 

The  slurry  consolidometer  (Figure  8.1)  consists  of  two  pieces  of  PVC  tube  with  an 
inside  diameter  of  524  mm  which  are  tighed  together  by  steel  rods.  The  entire  unit 
rests  on  a  frame  with  casters  attached  to  the  bottom  to  provide  mobility. 

The  uper  PVC  pipe  (see  Figure  8.1)  houses  a  piston  and  provides  additional  room 
to  hold  the  slurry  during  the  early  stage  of  slurry  consolidation.  The  bottom  com¬ 
partment  was  split  into  two  sections  longitudinally  and  a  steel  frame  was  mounted 
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Figure  8.1:  Cross  sectional  view  of  the  slurry  consolidometer. 
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around  each  piece  in  order  to  clamp  them  together  and  ultimately  take  them  apart 
once  the  soil  specimen  was  ready.  A  rubber  gasket  was  glued  to  the  cut  section  of 
the  lower  compartment  to  seal  it.  Sand  paper  was  attached  to  the  inside  of  the  lower 
compartment  to  prevent  slippage  of  the  membrane  during  consolidation.  The  mem¬ 
brane  extended  out  of  the  lower  compartment  and  acted  as  a  seal  between  the  two 
compartments. 

The  piston  consists  of  an  aluminum  plate  with  porous  plastic  attached  to  it  and 
an  aluminum  rod  screwed  into  it  to  provide  some  stability.  The  piston  is  moved  down 
as  air  is  pumped  into  a  number  of  rubber  truck  inner  tubes  that  are  joined  together 
with  3.2  mm  nylon  tubing.  The  bottom  aluminum  plate  is  also  the  top  of  the  piston 
for  the  calibration  chamber,  and  along  with  the  slurry  consolidometer  piston  pro¬ 
vides  drainage  at  the  top  and  bottom  of  the  specimen  during  slurry  consolidation. 
The  consolidometer  was  designed  for  an  initial  slurry  height  of  1186  mm  and  a  final 
height  of  813  mm  with  a  diameter  of  524  mm. 


8.2  The  Calibration  Chamber 

The  entire  calibration  chamber  setup  along  with  the  pressure  control  system  and  the 
computerized  data  logging  and  control  is  shown  in  Figure  8.2.  The  chamber  (Figure 
8.3)  was  desgned  to  house  a  524  mm  diameter  soil  specimen  with  a  miaximum  height 
of  813  mm.  A  constant  stress  or  zero  strain  condition  can  be  applied  in  the  vertical  and 
the  horizontal  direction.  The  piston  applies  the  vertical  pressure  and  provides  bottom 
drainage.  The  guide  shaft  keeps  the  piston  straight  as  well  as  provides  an  opening 
for  the  transducer  wires  and  drainage  lines.  The  guide  shaft  runs  through  a  linear 
ball  bearing  that  was  mounted  underneath  the  base  plate.  The  bottom  plate  is  fitted 
with  porous  plastic,  drainage  lines,  13  piezometers  which  are  attached  to  pressure 
transducers,  and  a  pressure  transducer  mounted  flush  with  the  porous  plastic  that 
measures  <7„,  the  vertical  stress.  The  piezometers  were  placed  in  such  a  way  as  to 
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measure  pore  pressure  from  different  radii  away  from  the  penetrometer.  In  Figure  8.4, 
the  piezometer  locations  are  shown  with  their  respective  radii  away  from  the  center 
of  the  penetrometer.  The  piezometers  were  made  of  1.65  mm  diameter,  0.23  mm  wall 
stainless  steel  tubings.  They  extended  up  250  mm  from  the  bottom  plate.  The  tip 
was  filtered  with  porous  plastic.  The  stem  of  the  stainless  steel  tubing  was  attached 
to  a  pressure  transducer.  All  piezometers  were  installed  and  flushed  with  deaired 
water  before  the  slurry  consolidation.  Upon  insertion  of  the  penetrometer  from  the 
top  of  the  specimen,  the  piezometers  will  not  interfere  with  the  soil  response.  This  is 
an  important  advantage  of  having  the  piezometers  extend  up  from  the  bottom  plate. 

The  chamber  has  two  stainless  steel  drums  (910  mm  high,  6.5  mm  thick)  which 
form  the  "double  wall”.  This  is  required  if  the  specimen  is  to  be  consolidated  under 
K0  conditions.  When  the  pressure  in  both  the  cavities  formed  by  thetwo  walls  is  set 
equal,  the  inner  wall  becomes  rigid. 

The  chamber  system  was  designed  to  handle  1000  kPa  maximum  pressure.  The 
pressure  was  regulated  through  a  control  panel  by  using  Fairchild  model  1000  regula¬ 
tors  for  manual  control  and  Fairchild  model  T5700  I/P  and  E/P  transducer  regulators 
for  computer  control.  Gages  and  pressure  transducers  were  used  to  monitor  the  pres¬ 
sure. 


8.3  Model  Piezoblade 

The  model  piezoblade  was  machined  from  a  single  piece  of  stainless  steel  and  was  117 
mm  long,  42  mm  wide  and  7mm  thick  (Figure  8.5).  Four  cavities  were  incorporated 
to  allow  pore  pressure  measurements  at  different  parts  of  the  blade.  The  four  cavities 
were  filled  with  a  porous  plastic  and  were  connected  to  transducer  ports  outside  the 
chamber  by  means  of  stainless  steel  tubes  (1.65  mm  OD,  0.23  mm  wall  thickness). 
These  tubes  were  laid  in  grooves  cut  on  the  face  of  the  blade  and  through  a  12.7  mm 
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Figure  8.2:  Calibration  chamber  control  system. 
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Figure  8.3:  Cross  sectional  view  of  the  calibration  chamber. 
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OD  hollow  stainless  steel  rod  914  mm  long.  The  grooves  on  the  blade  were  filled  with 
an  epoxy  to  give  a  smooth  finish  to  the  blade.  The  metal  rod,  which  was  used  to 
push  the  blade,  was  also  attached  to  the  blade  with  epoxy. 


8.4  Model  Dilatometer  (DMT) 

The  model  Dilatometer  (Figure  8.6)  has  the  same  dimensions  as  the  model  piezoblade. 
It  is  exactly  1/2  scale  of  the  Marchetti  flat  Dilatometer.  The  tip  apex  angle  is  20° 
and  it  has  a  30  mm  diameter  expandable  diaphram  which  was  instrumented  with  a 
strain  gage.  The  maximum  membrane  expansion  allowed  is  0.5  mm.  The  membrane 
itself  was  made  of  silicone  rubber. 

The  same  type  of  hollow  rod  used  for  pushing  the  piezoblade  was  attached  to  the 
model  Dilatometer.  A  stainless  steel  tube  used  for  pressure  flow  into  the  cavity  along 
with  the  wires  for  the  strain  gage  were  placed  through  this  rod. 


8.5  Blade  Housing 

It  was  very  important  to  have  the  soil  contained  inside  the  chamber  during  the  consol¬ 
idation  and  consequent  blade  insertion.  In  order  to  keep  the  back  pressure  on  the  soil 
specimen,  a  housing  had  to  be  developed  to  hold  the  blade  during  the  consolidation 
phase  of  the  procedure.  This  housing  is  shown  in  Figure  8.7. 

This  unit  fit  into  the  top  plate  of  the  calibration  chamber  and  into  the  top  drainage 
plate  (Figure  8.3).  Drainage  ports  were  provided  to  flush  the  unit  and  to  apply  the 
back  pressure  to  the  exposed  part  of  the  soil.  A  teflon  piece  acted  as  the  buffer  be¬ 
tween  the  soil  and  the  blade.  This  teflon  piece  had  a  groove  cut  into  it  so  that  the 
blade  and  the  rod  could  fit  through.  This  unit  also  acted  as  a  guide  for  the  rod.  A 
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Figure  8.5:  Geometry  and  dimensions  of  the  model  piezoblade. 
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Figure  8.6:  Geometry  and  dimensions  of  the  model  dilatometer. 
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Figure  8.7:  Cross  sectional  view  of  the  blade  housing. 


linear  ball  bearing  was  used  to  keep  the  rod  straight  during  the  penetration. 


8.6  Data  Logging  and  Electrical  Control  System 

A  Zenith  248  personal  computer  was  used  along  with  a  Meirabyte  DAC-02  board  and 
a  Keithley  Series  500  data  logging  system  to  read  the  pressure  transducers  and  to 
control  the  cell  and  piston  pressure  during  anisotropic  consolidation.  To  control  the 
penetration  of  the  model  penetrometers,  a  stepper  motor  driven  linear  actuator  was 
incorporated.  By  using  the  stepper  motor,  the  speed  and  exact  position  of  the  blade 
could  be  controlled. 

Hewlett-Packard  DC  power  supplies  were  used  to  supply  the  excitation  voltage  to 
all  the  transducers  and  the  strain  gages  in  the  model  Dilatometer.  An  uninterrupt- 
able  power  supply  system  was  employed  in  case  of  a  power  failure. 


8.7  Slurry  Preparation  and  Consolidation 

K50  clay  was  the  only  soil  used  in  all  the  calibration  chamber  tests.  The  purpose  of 
a  slurry  consolidometer  is  to  produce  uniform  and  reproducible  research  specimens  of 
cohesive  soils.  It  is  desireable  to  obtain  a  high  degree  of  homogeneity  both  physically 
and  chemically  and  also  to  produce  a  soil  specimen  with  a  known  stress  history. 

The  water  content  that  was  selected  for  the  slurry  was  chosen  by  considering  a 
number  of  factors  as  revealed  by  Krizek  and  Sheeran  (1970).  A  high  water  content 
will  make  the  slurry  easy  to  deair  and  easy  to  place  in  the  consolidometer.  However, 
a  low  water  content  minimizes  the  segregation  of  soil  particles  and  makes  the  consol¬ 
idation  time  shorter.  They  recommended  for  a  kaolinite  mixture  a  water  content  of 
1.5  to  2  times  the  liquid  limit,  where  the  actual  amount  is  found  by  trial  and  error. 
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For  the  mixture  used  for  this  research  work,  it  was  found  that  2  times  the  liquid  limit 
was  sufficient  to  satisfy  the  criteria  mentioned  above.  This  corresponds  to  an  initial 
water  content  of  60%. 

Deionized  and  deaired  water  was  used  for  the  slurry  mixture.  The  kaolinite  and 
sand  were  alternately  placed  in  the  water  a  little  at  a  time  to  provide  some  mixture  in 
the  beginning.  The  slurry  was  prepared  in  a  50  gallon  tank,  and  a  mixer  was  used  to 
make  the  slurry  physically  homogeneous.  The  slurry  was  scooped  out  of  the  mixing 
tank  and  poured  into  the  consolidometer.  By  pouring  all  the  slurry  in  the  same  place 
and  letting  it  flow  to  the  other  edges  the  air  was  forced  out.  The  air  that  was  added 
by  the  mixer  or  was  still  in  the  slurry  for  other  reasons  appeared  to  come  out  as  the 
slurry  was  slowly  poured  into  the  consolidometer.  This  method  was  also  used  when 
preparing  smaller  specimens,  which  were  used  for  the  stepped  loading  consolidation 
tests  and  the  triaxial  tests. 

The  slurry  was  consolidated  under  a  vertical  stress  of  138  kPa  in  the  slurry  con¬ 
solidometer.  It  took  approximately  30  days  to  consolidate  one  specimen.  At  the 
end  of  slurry  consolidation,  a  524  mm  diameter  clay  specimen  with  a  height  of  ap¬ 
proximately  800  mm,  encased  in  a  rubber  membrane  was  created.  The  specimen  was 
directly  transported  to  the  chamber  for  final  specimen  preparation  without  trimming. 


8.8  Calibration  Chamber  Consolidation  and  Pen¬ 
etration  Tests 

The  specimens  had  water  contents  of  27-28%  after  consolidation  in  the  slurry  con¬ 
solidometer.  The  Modified  Cam  Clay  model  used  as  part  of  the  strain  path  analyses 
assumes  isotropic  soil  behavior.  The  isotropic  consolidation  of  the  chamber  specimen 
therefore  provides  a  more  realistic  comparison  with  the  strain  path  studies.  However, 
natural  soil  deposit  is  generally  consolidated  anisotropically  by  the  overburden  pres- 
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sure.  The  idea  of  performing  chamber  tests  in  both  isotropically  and  anisotropically 
consolidated  specimens  was  to  first  validate  the  analytical  studies  and  then  evaluate 
the  effects  of  stress  anisotropy.  All  chamber  specimens  were  normally  consolidated. 

All  the  drainage  lines  were  flushed  prior  to  the  chamber  consolidation.  A  back 
pressure  was  used  during  consolidation.  For  isotropic  consolidation,  an  effective  con¬ 
solidation  stress  of  218  kPa  w’as  applied  to  the  cell  and  piston  while  maintaining 
drainage  from  the  top  and  bottom  of  the  specimen. 

For  anisotropic  consolidation,  a  K  value  of  0.52  was  utilized  as  in  the  case 

of  CAU  triaxial  tests.  A  procedure  similar  to  that  of  controlled  gradient  consolida¬ 
tion  tests  was  developed  for  the  anisotropic  chamber  consolidation.  The  piezometers 
installed  inside  of  the  specimen  had  their  tips  at  approximately  the  same  level  as 
the  center  of  the  specimen.  They  were  used  to  monitor  the  pore  pressure  at  that 
level  of  the  specimen  during  consolidation.  To  maintain  a  constant  gradient,  the  cell 
pressure  (cr/,)  was  adjusted  to  keep  the  pore  pressure  readings  from  these  piezometers 
at  69  kPa  above  the  back  pressure.  The  piston  pressure  (<r„)  was  adjusted  accord¬ 
ingly  to  maintain  a  K  value  of  0.52.  Drainage  against  the  back  pressure  was  allowed 
at  the  top  and  bottom  of  the  specimen.  The  double  drainage  enables  the  consoli¬ 
dation  to  be  completed  within  a  reasonable  period  of  time  while  a  gradient  of  18.5 
was  maintained  between  the  center  and  ends  of  the  specimen.  This  consolidation  can 
only  be  conducted  with  a  computer  servo-controlled  system.  The  consolidation  took 
aproximately  one  month  regardless  of  the  stress  anisotropy.  The  experimental  set  up 
allowed  one  specimen  being  consolidated  in  the  slurry  consolidometer  while  another 
one  was  in  the  chamber.  Overall,  a  test  was  conducted  every  45  days.  The  status  of 
consolidation  was  observed  by  monitoring  the  water  level  in  the  back  pressure  tank 
along  with  the  piezometer  or  pore  pressure  readings.  Representative  results  from 
these  measurements  are  shown  in  Figures  8.8  and  8.9. 

There  were  a  total  of  three  model  piezoblade  tests  performed,  the  penetration  rate 
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Tube  Reading, 
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Figure  8.8:  Consolidation  curve  found  by  tank  reading. 
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Figure  8.9:  Consolidation  curve  found  by  piezometer  reading. 
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in  the  first  piezoblade  test  was  faster  than  for  the  remainding  tests.  This  provided 
a  comparison  between  penetration  rates.  For  the  remainding  tests  the  blades  were 
pushed  at  a  constant  rate  of  2  cm/s ec. 

Three  model  Dilatometer  tests  were  performed.  Prior  to  the  insertion  of  the  model 
Dilatometer,  the  diaphram  was  inflated  in  the  blade  housing  against  the  back  pres¬ 
sure  and  then  deflated.  This  procedure  provided  the  measurement  of  the  membrane 
stiffness  needed  for  the  interpretation  of  the  test  result. 

During  the  first  minute  after  the  DMT  was  pushed  into  the  specimen,  a  full 
Dilatometer  inflation  test  was  performed  with  maximum  membrane  inflation  of  1.0 
mm.  Following  the  initial  full  expansion  tests,  cycles  of  inflating  the  membrane  to  just 
pass  lift  off  and  back  again  were  performed.  The  lift-off  pressure  is  generally  referred 
to  as  the  A  reading  and  the  pressure  where  the  membrane  returns  to  its  deflated 
position  is  called  the  C  reading  (Lutenegger,  1988).  It  is  important  to  note  here  that 
the  full  inflation  was  only  performed  initially  after  penetration.  Pore  pressure  was 
monitored  and  A  and  C  readings  were  taken  until  excess  pore  pressure  induced  by 
the  penetration  had  fully  dissipated. 


8.9  Experiments  After  the  Penetration  Test 

8.9.1  Moisture  Content  Determinations 

After  the  test  was  completed  moisture  content  samples  were  taken  from  the  specimen. 
The  samples  were  taken  at  different  levels  through  the  height  of  the  specimen  and 
at  J:.ferent  distances  away  from  the  path  of  the  blade.  Twenty  samples  each  were 
taken  roughly  25  mm,  150  mm,  305  mm,  457  mm,  and  610  mm  from  the  top  of  the 
specimen.  The  samples  were  taken  parallel  and  perpendicular  to  the  blade  path  on 
the  horizontal  plane  starting  25  mm  away  from  the  center  and  50  mm  away  from  the 
previous  sample  after  that.  In  general,  more  than  100  specimens  were  taken  for  water 
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content  tests  after  each  chamber  test. 


Upon  determination  of  the  moisture  contents  graphs  were  made  to  examine  the 
affect  of  blade  penetration  on  the  moisture  content.  It  was  expected  that  as  the  sam¬ 
ples  were  taken  closer  to  the  path  of  the  blade  the  moisture  content  would  decrease. 
The  results  are  summarized  and  discussed  in  the  next  chapter. 


8.9.2  Controlled  Gradient  Consolidation  Tests 

Along  with  moisture  content  tests,  soil  specimens  were  obtained  by  pushing  stainless 
steel  rings  (62  mm  ID,  17  mm  ht.)  into  the  chamber  specimen  for  controlled  gradient 
consolidation  tests.  These  specimens  were  taken  in  the  horizontal  direction  to  obtain 
the  corresponding  consolidation  characteristics.  The  same  controlled  gradient  consol¬ 
idation  testing  procedure  as  previously  described  was  used. 


8.10  Calibration  Chamber  Testing  Program 

A  total  of  six  flat-plate  penetration  tests  have  been  conducted  in  the  calibration 
chamber.  Table  8.1  furnishes  a  summary  of  the  piezoblade  (PIEZ1  to  PIEZ3)  and 
model  dilatometer  tests  (DMT1  to  DMT3).  The  folowing  sections  provide  a  detailed 
description  of  the  testing  procedure. 


164 


Test  No. 

Consolidation 

Back 

Pres.,  kPa 

Vert. 

Stress,  kPa 

Horiz. 

Stress,  kPa 

Penetration 

Rate,  cm/s 

PIEZl 

Iso. 

690 

896 

896 

6.0 

PIEZ2 

Iso. 

690 

896 

896 

2.0 

PIEZ3 

Aniso. 

690 

896 

797 

2.0 

DMT1 

Iso. 

207 

414 

414 

2.0 

DMT2 

Iso. 

552 

759 

759 

2.0 

DMT3 

Aniso. 

552 

759 

659 

2.0 

Table  8.1:  Summary  of  the  calibration  chamber  tests. 
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Chapter  9 


Calibration  Chamber  Test  Results 


Results  of  the  chamber  calibration  tests  are  presented  in  this  chapter.  These  include 
the  penetration  induced  pore  pressure,  the  pore  pressure  dissipation  after  interruption 
for  the  piezoblade  tests  and  the  expansion  curves  of  the  model  Dilatometer  tests.  The 
same  coordinate  system  as  used  in  presenting  the  numerical  data  (see  Figure  3.2)  is 
adopted.  The  half  thickness  of  the  blade  (R  =  3.5  mm)  was  used  to  normalize  linear 
dimensions.  When  appropriate,  stresses  and  pore  pressures  are  normalized  with  re¬ 
spect  to  vertical  consolidation  stress. 


9.1  Penetration  Induced  Pore  Pressure 

The  four  piezo  elements  (sensors)  mounted  on  the  piezoblade  enabled  excess  pore 
pressures  to  be  measured  at  different  locations  on  the  penetrometer  during  a  steady 
penetration.  In  test  PIEZl  (see  Table  8.1),  the  rate  of  penetration  was  6  cm/sec. 
At  this  penetration  rate,  large  negative  pore  pressure  readings  which  exceeded  the 
lower  limit  set  up  in  the  data  logging  system  was  registered.  Such  a  phenomenum 
has  not  been  reported  for  cone  penetration  tests  at  a  comparable  penetration  rate. 
The  complete  time  history  of  the  normalized  excess  pore  pressure,  A u/Pq  since  the 
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end  of  penetration  is  shown  in  Figure  9.1.  Despite  the  significant  difference  in  pen¬ 
etration  rate  and  excess  pore  pressure  developed  during  penetration,  the  maximum 
pore  pressure  recorded  in  PIEZl  is  very  close  to  that  in  PIEZ2  where  both  tests  were 
conducted  in  an  isotropically  and  normally  consolidated  specimen.  Figure  9.2  shows 
the  profiles  of  A u/Pq  developed  during  the  penetration  in  tests  PIEZ2  and  PIEZ3. 

A u  =  ur  -Uq  (9.1) 


where  ut  =  total  recorded  pore  pressure  and  u0  =  back  pressure.  A  penetration  rate 
of  2  cm/sec  was  used  in  these  two  tests  and  a  relatively  stable  excess  pore  pressure 
was  developed  at  150  mm  after  the  penetration  was  initiated.  In  test  PIEZ2  the 
clay  specimen  was  normally  consolidated  under  an  isotropic  condition,  in  test  PIEZ3 
the  consolidation  was  conducted  anisotropically  (see  Table  8.1).  For  both  tests,  the 
excess  pore  pressure  at  the  mid  height  of  the  blade  (sensors  3  and  4)  was  higher  than 
that  recorded  torward  the  end  of  the  blade  (sensors  1  and  2).  Also,  the  pore  pressure 
at  the  center  of  the  blade  (sensors  1  and  3)  was  higher  than  that  recorded  close  to 
the  edge  (sensors  2  and  4)  of  the  blade.  Sensor  No.  2  of  test  PIEZ3  was  accidentally 
disconnected  during  penetration  and  therefore  results  are  not  shown.  The  normal¬ 
ized  excess  pore  pressure,  A  u/Pq  from  all  the  sensors  ranged  from  1.5  to  2  during  the 
steady  penetration  in  test  PIEZ2.  The  correspoding  values  in  test  PIEZ3  varied  from 
0.9  to  1.4.  For  an  su/Pq  of  0.28  (see  Table  7.3)  as  indicated  in  the  reference  triaxial 
tests,  these  values  correspond  to  A u/su  ranging  from  5.4  to  7.1  for  test  PIEZ2  and 
from  3.2  to  5  for  test  PIEZ3.  These  values  are  well  within  the  range  of  those  predicted 
according  to  the  cavity  expansion  theory  (Vesic,  1972)  where 
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Figure  9.1:  Excess  pore  pressure  during  and  after  penetration  (PIEZ1). 
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The  penetration  induced  pore  pressure  in  the  isotropically  consolidated  specimen 
(PIEZ2)  was  higher  than  that  in  an  anisotropically  consolidated  specimen  (PIEZ3) 
by  as  much  as  43%.  According  to  Randolph  and  Wroth  (1979),  for  cylindrical  pen¬ 
etrometers 

Au  =  su  [ln{G/su)~  2 ln(z/R)}  (9.4) 

The  maximum  Au  is  at  the  boundary  of  the  penetrometer  where  z  =  R.  The  reference 
triaxial  tests  indicate  that 

=  2.5  (9 

for  specimens  normally  consolidated  to  218  kPa  and  sheared  under  axial  extension 
mode.  The  same  ratio  from  axial  compression  tests  was  0.83.  Soil  behavior  under 
axial  extension  is  similar  to  that  of  lateral  compression.  The  data  would  indicate 
therefore  that  the  soil  is  mostly  experiencing  a  lateral  compression  during  the  later 
stage  of  a  flat-plate  penetration  and  that  the  rigidity  index  in  the  lateral  direction  is 
the  predominant  factor  in  the  development  of  excess  pore  pressure. 

Figure  9.3  shows  a  comparison  of  A u/Pq  from  chamber  tests  and  those  from  SPM 
analyses  for  DMT  penetration.  It  should  be  emphasized  that  the  piezoblade  has  ex¬ 
actly  the  same  dimensions  as  the  model  DMT.  The  analytical  results  also  include 
A  u/Pq  determined  during  a  20°  cone  penetration.  Since  pore  pressure  can  only  be 
analytically  computed  along  the  x  =  0  plane,  the  chamber  data  shown  in  Figure  9.3 
represent  those  from  sensors  1  and  3.  The  figure  shows  that  the  analytical  values 
are  approximately  50%  of  those  obtained  in  test  PIEZ2  and  70%  of  test  PIEZ3.  The 
significant  difference  is  most  likely  due  to  the  way  MCC  model  handles  soil  unload 
and  reload.  As  indicated  in  Chapter  3,  the  shear  octahedral  strain  undergoes  a  pro¬ 
cess  of  unload  and  reload  as  the  soil  element  passes  the  tip  of  the  DMT  along  x  =  0 
plane.  The  initial  loading  and  unloading  of  the  soil  was  mainly  in  a  mode  of  triaxial 
compression  as  predicted  in  SPM  studies.  The  MCC  model  considers  unload-reload 


( suln(G/su))CTV 
{suln(G/su))CAU 


as  elastic  and  therefore  there  is  no  shear  stress  induced  pore  pressure  during  this 
stage.  This  is  clearly  not  the  case  as  indicated  in  reference  triaxial  tests  with  unload 
reload  cycles.  The  data  base  regarding  the  post  failure  unload-reload  behavior  is  very 
limited.  It  is  clear  that  extensive  experimental  studies  are  needed  to  establish  such  a 
data  base  before  the  analytical  model  can  be  improved. 

Despite  the  discrepancies,  the  analytical  study  does  correctly  predict  the  trend  of 
excess  pore  distribution  which  shows  that  A u  decreases  towards  the  end  of  the  blade. 
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Figure  9.2:  Profiles  of  excess  pore  pressure  during  piezo  blade  penetration. 
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Figure  9.3:  Comparison  with  analytical  results. 
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9.2  Pore  Pressure  Distribution  and  Dissipation 


The  following  discussion  will  concentrate  on  results  from  tests  PIEZ2  and  PIEZ3 
because  of  their  similarity  in  penetration  rate  and  a  desirable  back  pressure  applied 
during  the  test.  Piezometers  were  strategically  located  within  the  clay  specimen  to  al¬ 
low  the  pore  pressure  distribution  to  be  monitored.  Since  the  piezometers  were  made 
of  relatively  flexible  stainless  steel  needles,  they  could  be  moved  by  the  surrounding 
soil  during  the  specimen  consolidation.  The  exact  location  of  each  piezometers  were 
therefore  measured  by  carefully  exposing  them  after  the  test  was  completed.  The 
"as-measured”  locations  of  all  the  piezometers  along  with  those  of  the  sensors  (No.  3 
and  4)  on  the  piezoblade  are  shown  in  Figures  9.4  and  9.5.  It  should  be  noticed  that 
the  tip  of  the  piezometers  was  at  the  same  level  as  sensors  3  and  4  on  the  piezoblade. 
Although  the  piezometers  were  carefully  flushed  with  deaired  water  and  the  specimen 
was  under  a  back  pressure  of  690  kPa  during  the  test,  an  apparent  delayed  response 
was  noticed  in  all  the  piezometers.  This  could  be  caused  by  a  combined  effect  of  soil 
compressibility  and  the  compliance  of  the  piezometer /pressure  transducer  system.  In 
any  case,  stable  readings  were  obtained  typically  at  10  seconds  after  the  interruption 
of  the  penetration.  Tables  9.1  and  9.2  show  all  the  piezometer  readings  in  terms  of 
normalized  excess  pore  pressure  (A u/Pq)  at  selected  time  intervals  after  the  pene¬ 
tration  ended.  To  visualize  the  distribution  of  the  pore  pressure,  contour  lines  were 
created  by  interpolating  these  piezometer  readings  using  a  commercial  software  called 
SURFER.  Contours  resulted  from  these  interpolations  (Figures  9.6  to  9.17)  represent 
the  distribution  of  excess  pore  pressure  in  the  xz  plane  and  at  the  same  level  as  pore 
pressure  sensors  3  and  4  (y/R  =  13).  Based  on  the  shapes  of  the  contours,  it  appears 
that  the  pore  pressure  was  fairly  uniform  across  the  face  of  the  piezoblade  (in  the 
x  direction)  during  dissipation  for  test  PIEZ2  (isotropically  consolidated  specimen). 
For  test  PIEZ3  (anisotropicaly  consolidated  specimen)  the  pore  pressure  dissipation 
follows  an  axisymmetric  pattern. 
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Figure  9.4:  Position  of  piezometers  after  test  (PIEZ2). 
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Figure  9.5:  Position  of  piezometers  after  test  (PIEZ3). 
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Time 

Piezometer  No. 

sec. 

1 

2 

3 

D 

5 

D 

D 

8 

9 

11 

12 

13 

10 

0.04 

0.06 

0.06 

0.04 

0.02 

0.01 

0.01 

0.08 

0.52 

0.02 

0.01 

0.01 

50 

0.10 

0.12 

0.12 

0.08 

0.03 

0.01 

0.01 

0.25 

0.86 

0.05 

0.02 

0.01 

100 

0.16 

0.19 

0.49 

0.13 

0.05 

0.00 

0.01 

0.34 

0.97 

0.09 

0.02 

0.02 

1000 

0.51 

0.54 

0.49 

0.37 

0.20 

0.08 

0.08 

0.72 

0.69 

0.31 

0.10 

0.13 

3000 

0.47 

0.42 

0.38 

0.32 

0.26 

0.21 

0.20 

0.45 

0.43 

0.29 

0.21 

0.23 

6000 

0.34 

0.32 

0.31 

0.29 

0.20 

0.26 

0.26 

0.33 

0.33 

0.27 

0.26 

0.26 

Table  9.1:  Table  of  piezometer  readings  for  PIEZ2(Au/P'). 


Time 

Piezometer  No. 

sec. 

1 

2 

3 

4 

5 

6 

fl 

8 

9 

10 

11 

12 

13 

10 

0.23 

0.08 

0.03 

0.01 

0.00 

0.00 

0.00 

0.13 

0.04 

0.25 

0.00 

0.00 

0.00 

50 

0.71 

0.44 

0.16 

0.03 

0.00 

0.00 

0.00 

0.65 

0.17 

0.70 

0.04 

0.01 

0.00 

100 

0.82 

0.66 

0.27 

0.06 

0.00 

0.00 

0.00 

0.81 

0.24 

0.79 

0.06 

0.01 

0.01 

1000 

0.59 

0.55 

0.44 

0.21 

0.00 

0.02 

0.03 

0.53 

0.18 

0.45 

0.19 

0.07 

0.03 

3000 

0.29 

0.27 

0.24 

0.17 

0.03 

0.06 

0.07 

0.27 

0.11 

0.24 

0.15 

0.10 

0.07 

6000 

0.17 

0.17 

0.16 

0.13 

0.06 

0.09 

0.09 

0.17 

0.09 

0.17 

0.12 

0.10 

0.09 

Table  9.2:  Table  of  piezometer  readings  for  PIEZ3(A uj P'0). 
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Contour  of  pressure  6000  sec.  after  interruption  (PIEZ2). 
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Figure  9.14:  Contour  of  pressure  100  sec.  after  interruption  (PIEZ3). 
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9.2.1  Penetration  Induced  Excess  Pore  Pressure 


Figure  9.18  shows  the  excess  pore  pressure  in  the  z  direction  on  the  x  =  0  plane.  The 
figure  was  established  based  on  the  values  depicted  in  the  contours  at  10  seconds  after 
the  penetration  was  interrupted.  The  x  =  0  plane  was  selected  as  this  is  typically 
where  the  excess  pore  pressure  is  monitored  in  field  flat-plate  penetration  tests.  The 
amount  of  consolidation  in  10  seconds  is  minimal  and  therefore  the  distribution  shown 
in  Figure  9.18  is  considered  the  same  as  immediately  after  the  penetration  or  at  the 
beginning  of  pore  pressure  dissipation.  According  to  this  figure,  the  penetration  in¬ 
duced  excess  pore  pressure  diminishes  within  a  distance  of  20R  (R  =  half  thickness 
of  the  piezoblade).  This  is  significantly  less  than  that  predicted  using  the  equation 
proposed  for  axisymmetric  penetrations  (Randolph  and  Wroth,  1979).  In  that  sense 
the  pore  pressure  measurements  in  the  chamber  agree  well  with  the  SPM  analyses. 
The  flat-plate  penetration  induced  shear  octahedral  strain  diminishes  within  a  short 
distance  (as  opposed  to  cone  penetration). 

In  an  attempt  to  improve  the  analytical  description  of  the  pore  pressure  distri¬ 
bution,  the  concept  proposed  by  Liam  Finn  (1963)  was  adopted.  Considering  the 
flat-plate  penetration  as  a  boundary  value  problem  of  a  rigid  plate  translating  later¬ 
ally  in  an  elastic  half  space,  the  normal  stresses  along  the  center  line  (x  =  0)  of  the 
rigid  plate  are 


b  2  by2 

b2  +  z3  +  (b2  +  y2)2 

d  E  b  2  by2 

bn  (1  -  1/2)  [&2  +  “  ( b 2  -I-  2/2)2 


d  E 
bn  (1  —  I/2) 


where 

d  =  R,  the  half  thickness  of  the  flat-plate,  and 
b  =  half  width  of  the  flat  plate 


(9.6) 

(9.7) 
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The  excess  pore  pressures  generated  may  be  estimated  by  assuming  that,  under 
undrained  conditions,  the  mean  effective  stress  remains  constant.  Thus,  the  excess 
pore  pressure  is  equal  to  the  change  in  mean  total  stress,  giving 


A  1/A  A  a  \  1  /  A  a  \  d  E  [6 

Au  =  j(A<r,  +  Aav  +  At?*)  =  -(Ac,  +  A<t2)  = 


(9.8) 


The  ratio  of  A u  at  any  given  distance  z  away  from  the  face  of  the  blade  over  that  on 
the  face  of  the  blade,  A u0,  is  therefore 

Au  (b7  +  z7) 


Au0  (b*  +  R7Y 


(9.9) 


By  setting  A u0  equal  to  that  measured  from  sensor  3  at  the  end  of  penetration,  the 
initial  excess  pore  pressure  distribution  can  be  estimated  based  on  EQ.  9.9.  Accord¬ 
ing  to  EQ.  9.9  the  pore  pressure  distribution  is  independent  of  the  rigidity  of  the  soil 
provided  it  is  elastic  and  A u0  is  given.  In  field  practice,  A ua  is  the  measured  value 
from  the  piezo  unit  immediately  after  the  interruption  of  the  penetration.  Figure 
9.19  plots  the  predicted  pore  pressure  distribution  for  both  the  isotropically  (PIEZ2) 
and  anisotropically  (PIEZ3)  consolidated  specimens.  In  general,  the  predicted  exces¬ 
sive  pore  pressure  is  significantly  higher  than  those  according  to  chamber  piezometer 
readings.  Obviouly,  it  could  be  argued  that  soil  is  not  elastic  and  that  a  plastic  region 
is  typically  formed  around  the  penetrometer  which  significantly  reduces  the  area  of 
pore  pressure  increase  (Vesic,  1972).  Also,  EQ.  9.9  assumes  plane  strain  expansion 
whereas  the  pore  pressures  below  and  above  the  flat-plate  are  much  lower  in  real 
tests  and  therefore  allows  pore  pressure  to  dissipate  in  the  vertical  direction.  Despite 
the  above  problems,  EQ.  9.9  does  offer  a  basis  of  describing  the  initial  pore  pressure 
distribution  and  has  a  physical  meaning.  The  equation  shows  that  when  z  is  small 
(relative  to  b  value),  Au  changes  slowly  with  z  as  in  the  case  of  large  loading  area 
over  a  thin  soil  layer.  As  z  increases,  the  effect  of  b  diminishes  and  the  pore  pressure 
distribution  essentially  follows  the  ratio  of  ( z/R )2.  Following  this  concept  EQ.  9.9 
can  be  modified  as 

Au  _  (frm  +  zn) 

A u0  ~  (6m  +  Rn)' 
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where  m  and  n  are  variables  to  be  determined  from  curve  fitting  the  chamber  data. 
The  values  of  m  and  n  should  reflect  the  effects  of  the  three-dimensional  and  non¬ 
elastic  nature  of  the  flat-plate  penetration.  Also,  EQ.  9.10  forms  the  basis  of  a 
semi-empirical  approach  for  the  interpretation  of  pore  pressure  dissipation  data  to  be 
presented  later.  Results  of  the  curve  fitting  show  that  for  the  isotropically  consoli¬ 
dated  specimen  (PIEZ2) 


Au  (618  -f  z2-6) 
Au0  ~  (61-8  +  R2-6) 

and  for  anisotropically  consolidated  specimen  (PIEZ3) 

Au  ( b 27  +  z30) 
Au0  ~  (ft2-7  +  R3  0) ' 


(9.11) 


(9.12) 


Curves  based  on  Eqs.  9.11  and  9.12  are  also  plotted  in  Figure  9.18. 


9.2.2  Pore  Pressure  Dissipation  and  Determination  of  Con¬ 
solidation  Characteristics 

Figures  9.21  and  9.22  show  the  dissipation  of  excess  pore  pressure  after  the  penetration 
ended.  The  rate  of  pore  pressure  dissipation  is  governed  mostly  by  the  coefficient  of 
consolidation  of  the  surrounding  soil.  Attempts  have  been  made  to  determine  the  co¬ 
efficient  of  consolidation  based  on  the  pore  pressure  dissipation  record  from  flat-plate 
penetration  tests  (Kabir  and  Lutenegger,  1990).  All  these  methods  assume  flat-plate 
penetration  as  axisymmetric  and  procedures  originally  developed  for  cylindrical  pen¬ 
etrometers  have  been  followed.  The  success  of  such  methods  has  been  limited.  One 
important  drawback  of  the  current  approach  is  the  lack  of  knowledge  regarding  the 
distribution  of  pore  pressure  around  the  flat-plate  at  different  stages  of  pore  pressure 
dissipation. 


Analytical  studies  presented  previously  have  indicated  that  a  flat-plate  penetra¬ 
tion  is  a  three-dimensional  process.  Ideally,  a  three-dimensional  solution  should  be 
provided  to  properly  interpret  the  pore  pressure  dissipation  data  and  to  determine  the 
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coefficient  of  consolidation.  However,  the  coefficient  of  consolidation  is  influenced  by 
the  state  of  stress  and  could  vary  by  as  much  as  one  order  of  magnitude  as  indicated  in 
the  reference  consolidation  tests.  In  addition,  a  three-dimensional  distribution  of  ini¬ 
tial  excess  pore  pressure  would  be  necessary  for  such  a  solution.  It  would  be  extremely 
difficult  to  obtain  such  information,  even  in  a  calibration  chamber.  Therefore,  some 
simplifications  can  be  justified  to  determine  the  "approximate”  coefficient  of  consol¬ 
idation  in  a  reasonable  and  practical  fashion.  An  evaluation  of  the  pore  pressure 
contours  shown  in  Figures  9.6  to  9.17  shows  that  for  the  isotropically  consolidated 
specimen  (i.e.,  in  test  PIEZ2),  pore  pressure  was  relatively  uniform  along  the  face 
of  the  plate  throughout  the  dissipation  process.  For  the  anisotropically  consolidated 
specimen  (PIEZ3)  the  pore  pressure  distribution  follows  approximately  an  axisym- 
metric  pattern  during  dissipation.  If  consolidation  is  two-dimensional  (independent 
of  y  coordinate)  and  pore  pressure  is  monitored  at  the  center  of  the  flat-plate  (x=0, 
z=R),  the  pore  pressure  dissipation  should  be  governed  by  a  consolidation  process 
which  is  between  one-dimensional  (in  z  direction  only)  and  axisymmetric.  These  are 
further  simplified  as  one-dimensional  problems  with  a  governing  equation  as 


d2u  du 

Chd?  =  ~dt 
Ch.  d  ( rdu\ 
r  dr  \  dr  ) 


( one  dimensional ) 


du 

di 


( axisymmetric ) 


where 

Ch  =  horizontal  coefficient  of  consolidation 
u  =  excess  pore  pressure 
t  =  time  elapsed 


(9.13) 

(9.14) 


To  further  explore  the  possibility,  a  finite  difference  program  was  developed  to 
solve  the  above  equations.  Eqs.  9.11  and  9.12  were  used  to  determine  the  initial 
excess  pore  pressure  distribution.  Boundary  conditions  shown  in  Figure  9.23  were 
used  in  the  finite  difference  computation.  They  consider  the  penetrometer  surface 
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and  the  specimen  outer  boudary  as  impermeable  which  are  comparable  to  the  cases 
of  chamber  testing.  It  was  assumed  that  the  consolidation  was  linear  and  uncoupled 
so  that  c/i  and  the  total  stress  remain  constant.  Computational  results  in  terms  of 
A u/p'0  at  the  flat-plate  surface  (x=0,  z=R)  versus  time  factor  r(cht/R 2)  are  shown  in 
Figure  9.24.  Notice  that  A u  does  not  necessarily  reach  0  at  large  r.  This  is  because 
of  the  limited  and  impermeable  outer  boundary  used  in  the  analyses  which  is  also  the 
case  in  chamber  testing  (see  Figures  9.21  and  9.22). 


The  Ch  value  can  be  determined  by  comparing  the  pore  pressure  dissipation  data 
(i.e.,  Figures  9.21  and  9.22)  and  the  analytical  solution  (i.e.,  Figure  9.24).  t50  is  taken 
as  the  time  elapsed  when  50%  of  the  excess  pore  pressure  on  the  flat-plate  surface  has 
dissipated  according  to  the  dissipation  data.  The  corresponding  r50  is  selected  from 
the  analytical  solution.  Ch  is  then  calculated  as 


TsoR2 


t50 


(9.15) 


Table  9.3  shows  the  estimated  Ch  values  following  the  above  procedure.  In  general,  Ch 
values  based  on  axisymmetric  solutions  are  higher  than  those  from  one-dimensional 
solutions.  The  results  of  controlled  gradient  consolidation  tests  performed  on  samples 
taken  from  the  chamber  specimens  are  summarized  in  Table  9.4.  For  the  isotropically 
consolidated  specimen,  the  estimated  Ch  is  at  least  one  order  of  magnitude  smaller 
than  those  in  Tables  9.4  and  Cy  values  in  Table  7.2.  For  the  anisotropically  consoli¬ 
dated  specimen,  the  estimated  Ch  based  on  axisymmetric  solution  is  within  the  order 
of  magnitude  but  close  to  the  lower  bound  of  the  corresponding  values  in  Table  9.4. 
The  SPM  study  indicated  a  significant  decrease  of  pore  pressure  between  the  tip  of 
the  flat-plate  (y=0)  and  the  location  of  sensor  3  (y = 1 3)  which  is  typically  where  the 
pore  pressure  is  monitored  in  the  field  piezoblade  tests.  Because  of  this,  there  is  also 
a  tendency  for  the  pore  pressure  to  increase  towards  the  center  of  the  flat-plate  as  the 
pore  pressure  dissipates  away  from  it  after  an  interruption  of  penetration.  This  is  a 
likely  reason  as  to  why  typically  excess  pore  pressures  from  piezoblade  tests  dissipate 
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much  slower  than  those  obtained  from  piezocone  tests.  For  piezocone  tests,  the  pore 
pressure  is  generally  monitored  within  the  region  of  the  cone  tip  where  the  excess 
pore  pressure  is  usually  the  highest. 

It  appears  that  for  isotropically  consolidated  specimens,  a  reasonable  analysis  of 
pore  pressure  dissipation  around  a  flat-plate  would  have  to  consider  the  components 
in  the  y  direction,  ^or  anisotropically  consolidated  soil,  a  lower  bound  C/,  can  be 
estimated  considering  the  pore  pressure  dissipation  as  axisymmetric  and  the  initial 
excess  pore  pressure  distribution  follows  that  described  by  EQ.  9.12. 


Consolidation 

Ck,  10  3cm2/sec 

One-Dimensional 

Axisymmetric 

Isotropic 

0.48 

0.20 

Anisotropic 

0.50 

1.50 

Table  9.3:  Estimated  c k  from  the  piezoblade  tests. 


Tests 

C^,10  3cm2/sec 

Cc 

Reload 

Virgin  Loading 

Isotropic 

6.9  to  67.4 

4.6  to  8.6 

0.15  to  0.26 

Anisotropic 

11.1  to  72.5 

4.5  to  9.7 

0.19  to  0.23 

Table  9.4:  Summary  of  the  controlled  gradient  consolidation  tests  conducted  on  cham¬ 
ber  specimens. 
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Figure  9.20:  Dissipation  results  from  PIEZ2. 
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Figure  9.22:  Boundary  conditions  used  in  the  finite  difference  computation. 
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Figure  9.23:  Results  from  finite  difference  computation. 
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9.3  Results  of  the  Model  DMT  Tests 


Results  from  two  model  DMT  tests  are  presented  in  the  following  sections.  In  test 
DMT2,  the  diaphragm  was  ruptured  during  penetration  and  therefore  no  results  wil 
be  presented.  Test  DMT1  was  conducted  in  an  isotropically  consolidated  specimen 
and  DMT3  was  performed  in  an  anisotropically  consolidated  specimen.  In  both  tests 
the  diaphragm  was  expanded  shortly  after  the  DMT  penetration  was  interrupted.  In 
test  DMT1  the  diaphragm  was  inflated  to  roughly  0.5  mm  and  then  vented  to  its  orig¬ 
inal  position.  A  second  inflation  to  1  mm  was  then  performed  with  the  subsequent 
return  to  the  original  position.  These  were  the  only  full  inflation  tests  that  were  per¬ 
formed  in  DMTl.  In  test  DMT3,  the  diaphragm  was  inflated  to  roughly  0.5  mm  and 
then  vented  to  its  original  position,  no  further  full  expansion  was  applied  until  the 
excess  pore  pressure  was  fully  dissipated.  Notice  that  DMTl  was  performed  under 
a  back  pressure  of  218  kPa  and  a  back  pressure  of  550  kPa  was  applied  in  DMT3. 
Figures  9.25  and  9.26  show  the  applied  pressure  versus  diaphragm  expansion  from 
these  two  tests  and  stiffness  of  the  diaphragm  membrane  according  to  the  calibration 
tests  performed  in  the  blade  housing  unit  under  the  back  pressure.  For  each  DMT 
test,  a  set  of  P0,  Pi  and  P 2  pressures  were  computed  where 
P0  =  pressure  required  to  expand  the  diaphragm  off  the  face  of  the  DMT  blade 
P\  =  pressure  at  which  the  diaphragm  inflated  1  mm  away  from  the  face  of  the  blade 
P2  =  pressure  at  which  the  diaphragm  regains  its  contact  with  the  DMT  blade  fol¬ 
lowing  an  expansion. 

The  above  pressures  were  corrected  by  subtracting  the  membrane  stiffness  accord¬ 
ing  to  the  calibration  curves  which  included  back  pressure. 

While  the  pore  pressure  in  the  specimen  was  dissipating  a  succession  of  diaphragm 
inflations  and  deflations  were  performed  to  obtain  a  series  of  equivalent  DMT  P0  and 
P2  readings.  The  inflation  of  the  diaphragm  was  just  enough  to  find  the  lift-off  and 
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the  return  pressure.  Care  was  taken  not  to  induce  excessive  membrane  expansion. 


Because  of  the  different  dimensions  of  the  model  DMT  used  in  this  research  and 
a  conventional  DMT  blade,  the  interpretation  of  the  test  had  to  be  modified.  The 
value  of  Pi  was  considered  to  be  the  pressure  at  which  the  diaphragm  inflated  0.5 
mm  away  from  the  face  of  the  blade  and  not  1  mm.  This  number  was  chosen  because 
the  model  blade  dimensions  including  the  diaphragm  are  all  exactly  half  of  the  con¬ 
ventional  DMT  dimensions. 

Table  9.5  summarizes  the  parameters  derived  from  the  two  model  DMT  tests 
following  the  procedure  established  by  Marchetti  (1980).  According  to  Marchetti’s 
soil  classification  scheme  the  Id  values  from  both  tests  would  indicate  that  the  soil 
tested  is  clay.  This  result  correlates  well  to  the  actual  soil  type  which  was  silty-clay. 
Marchetti  (1980  )  suggested  that 

OCR  =  (Q.5KDy-5€  (9.16) 


Test  No. 

P0,  kPa 

Pi,  kPa 

P2,  kPa 

Id 

Kd 

Ep,  kPa 

UD 

DMT1 

396 

436 

285 

0.10 

1.91 

1532 

0.79 

DMT3 

309 

337 

226 

0.09 

1.49 

972 

0.73 

Table  9.5:  Summary  of  DMT  derived  Parameters. 

According  to  EQ.  9.16  the  OCR  value  is  0.93  for  DMT1  and  0.63  for  DMT3.  The 
actual  soil  was  normally  consolidated  with  an  OCR  of  1.  Considering  the  significant 
differences  in  dimensions  between  the  model  and  the  Marchetti  DMT,  these  empirical 
rules  perform  relatively  well. 


202 


Figures  9.27  and  9.28  show  a  comparison  between  the  P0  and  P2  readings  and 
the  excess  pore  pressures  taken  from  the  piezoblade  tests  (PIEZ2  and  PIEZ3).  All 
variables  are  normalized  with  respect  to  the  preconsolidation  pressure  P'a.  In  the  cor¬ 
responding  piezoblade  tests,  the  same  boundary  conditions  and  stress  history  were 
applied  as  in  the  model  DMT  tests.  The  same  excess  pore  pressures  measured  in  the 
piezoblade  tests  are  therefore  expected  to  occur  during  the  model  DMT  tests.  For 
tests  performed  both  in  the  isotropically  and  anisotropicaly  consolidated  specimens, 
Pq  was  higher  than  the  excess  pore  pressure  in  the  very  beginning  of  the  pore  pressure 
dissipation  process.  This  indicates  that  the  total  stress  around  the  DMT  is  somewhat 
higher  than  the  pore  pressure.  During  the  majority  of  the  pore  pressure  dissipation, 
Pq  remained  close  to  the  excess  pore  pressure  readings.  Towards  the  end  of  dissi¬ 
pation,  however,  Pq  was  higher  than  the  corresponding  excess  pore  pressures.  This 
further  confirms  the  existence  of  a  compressive  effective  stress  adjacent  to  the  DMT. 
The  P2  values  follow  closely  or  slightly  below  the  excess  pore  pressure  readings.  A 
suction  could  be  developed  during  deflation  of  the  DMT  diaphragm.  This  is  most 
likely  the  reason  for  the  lower  P2  readings  from  test  DMT3.  In  any  case,  the  trend  of 
P2  readings  versus  time  is  very  close  to  that  of  the  excess  pore  pressure  dissipation. 
Therefore,  it  would  be  possible  to  use  P2  readings  in  the  determination  of  q,  value. 

Figures.  9.29  to  9.32  show  the  moisture  content  distribution  down  the  front  and 
side  of  the  flat-plate  path  from  four  separate  specimens  all  isotropically  consolidated. 
There  is  no  marked  differences  in  the  moisture  contents  taken  from  the  four  specimens. 
This  would  indicate  the  reasonable  repeatability  of  specimen  preparation  procedure 
adopted  in  the  research  project. 

Figures.  9.33  and  9.34  show  the  moisture  content  distribution  from  the  front  and 
side  of  the  flat-plate  of  an  anisotropically  consolidated  specimen.  Ail  the  profiles 
indicated  a  consistent  decrease  of  moisture  content  by  as  much  as  2%  adjacent  to 
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the  DMT  blade.  This  corresponds  to  a  change  of  void  ratio  of  0.054.  If  flat-plate 
penetration  induces  strictly  a  one-dimensional  lateral  compression,  the  reduction  of 
such  void  ratio  would  mean  an  increase  of  lateral  consolidation  stress  on  the  order 
of  200  to  300  kPa  based  on  the  Cc  values  from  consolidation  tests  on  samples  taken 
from  the  chamber  specimens  (Table  9.4).  These  are  significantly  higher  than  the  P0 
readings  after  the  pore  pressure  dissipation.  It  is  possible  therefore  that  significant 
creeping  had  occurred  in  the  process  of  pore  pressure  dissipation. 
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Figure  9.24:  Expansion  curves  of  test  DMTl. 
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Figure  9.25:  Expansion  curves  of  test  DMT3. 
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Figure  9.26:  Normalized  P0  and  P2  readings  (DMT1)  and  piezoblade  dissipation 
(PIEZ2)  data. 
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Figure  9.27:  Normalized  Pq  and  Pi  readings  (DMT3)  and  piezoblade  dissipation 
(PIEZ3)  data. 
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Figure  9.30:  Moisture  contents  from  isotropically  consolidated  specimens  (down  the 
front  of  the  blade). 
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Figure  9.31:  Moisture  contents  from  isotropically  consolidated  specimens  (down  the 
side  of  the  blade). 
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Figure  9.32:  Moisture  contents  from  anisotropically  consolidated  specimens  (front  of 
the  blade). 
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Figure  9.33:  Moisture  contents  from  anisotropically  consolidated  specimens  (side  of 
the  blade). 
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Chapter  10 


Comparisons  with  Field  Tests 


A  comparison  was  made  of  the  results  obtained  from  the  various  flat-plate  penetrom¬ 
eters  by  conducting  field  tests  which  included  DMT,  AVStepped  Blade  and  tapered 
blade  tests.  The  test  site  is  located  east  of  the  town  of  Massena,  New  York,  which  is 
within  the  southern  extent  of  the  Champlain  sea.  The  soils  consist  of  approximately 
12  m  of  a  marine  clay  deposit  underlain  by  glacial  till.  Below  a  highly  weathered 
and  fissured  crust  of  1.5  m  thick,  the  marine  clay  is  soft  and  lightly  overconsoli¬ 
dated  (OCR  =  1.5  to  3)  .  In  addition  to  conventional  laboratory  classification  tests, 
reference  insitu  tests  performed  at  the  test  site  included  Self-Boring  Pressuremeter 
(SBPM),  Full-Displacement  Pressuremeter  (Pencel),  Field  Vane  Shear  (FVT)  and 
Piezocone  (CPTU)  tests.  Figure  10.1  shows  profiles  of  the  available  test  results. 

The  interpretation  of  the  /^-Stepped  and  tapered  blade  followed  the  procedures  by 
Handy,  et  al.  (1982)  which  uses  linear  regression  to  fit  the  data  to  EQ.  1.4  with  as¬ 
sumed  value  a  =  1.  A  coefficient  of  correlation  is  also  computed  during  the  linear 
regression.  The  results  are  shown  in  Figure  10.2  and  plots  of  the  normalized  pressure 
versus  blade  thickness  for  all  the  data  axe  included  in  Figure  10.3.  The  DMT  results 
in  terms  of  the  net  pressures  at  *he  beginning  (P0DMT)  and  end  of  the  membrane 
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expansion  (P^MT)  versus  depth  are  plotted  in  Figure  10.4  along  with  the  SBPM  and 
Pencel  lift  off  and  limit  pressures.  In  essence,  the  field  flat-plate  test  results  support 
the  SPM  studies  substantially.  The  pressure  readings  from  the  first  sensor  are  con¬ 
sistently  higher  than  those  from  the  second  one  in  most  of  the  stepped  and  tapered 
blade  tests.  These  results  are  similar  to  those  obtained  at  other  sites  (Lutenegger 
and  Timian,  1986a;  1986b)  using  the  same  AVStepped  Blade  and  are  consistent  w*ith 
results  obtained  using  a  widened  blade  (Handy,  et  al.  ,  1987).  These  penetrome¬ 
ters  have  relatively  large  apex  angles.  Strain  path  analyses  indicated  that  there  is 
a  large  reversal  of  the  ^-strain  as  the  soil  passes  the  blade  tip  Computations  of 
pore  pressure  and  total  stress  presented  earlier  indicated  that  there  is  sharp  increase 
and  subsequent  decrease  of  pore  pressure  at  the  tip  of  a  penetrometer,  regardless  of 
its  geometry.  Due  to  the  proximity  of  the  first  pressure  sensor  to  the  blade  tip,  the 
sensor  is  likely  to  be  influenced  by  such  a  high  pore  pressure  thus  yield  high  pres¬ 
sure  readings.  This  abnormality  is  responsible  for  some  of  the  negative  slopes  for  the 
AVStepped  Blade  tests  and  undesirable  coefficients  of  correlation  for  all  as  shown  in 
Figure  10.2.  Neglecting  the  first  pressure  reading  is  not  preferable  as  it  would  jeopar¬ 
dize  the  statistical  validity  of  the  data  especially  for  the  AVStepped  Blade.  The  test 
data  shown  in  Figures  refkoblade  and  10.3  indicate  that  for  the  tapered-blade,  the 
measured  lateral  pressures  are  rather  insensitive  to  the  blade  thickness  in  this  clay. 
This  is  consistent  with  the  analytical  studies  which  indicate  that  the  pore  pressure 
and  total  lateral  stress  (crzz)  remain  relatively  constant  beyond  the  tip  of  the  blade, 
as  the  blade  thickness  increases. 

The  flexible  membrane  on  the  DMT  extends  to  approximately  x  =  4.2 R  (30  mm 
radius).  Based  on  SPM  analyses,  the  jEVstrain  within  the  membrane  area  should 
be  between  those  of  plane  strain  and  CPT  penetrations.  The  PqMT  value  should 
therefore  reflect  an  effect  between  plane  strain  and  cylindrical  cavity  expansions. 
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This  is  consistent  with  the  PqMT  values  shown  in  Figure  10.4  which  are  very  close  to 
the  SBPM  limit  pressures.  Similar  results  were  also  reported  by  Lutenegger  (1988) 
for  other  DMT  tests  performed  in  the  same  marine  clay  deposit. 
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:  Soil  profiles  at  the  test  site. 
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Figure  10.3:  Normalized  stepped  and  tapered  blade  test  results. 


Figure  10.4:  SBPM  and  DMT  test  results. 
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Chapter  11 
Conclusions 


In  this  research,  a  numerical  technique  capable  of  computing  strain  paths  for  three- 
dimensional  penetrometers  was  developed.  It  provided  a  systematic  way  of  elucidat¬ 
ing  the  soil  strain  and  displacement  induced  by  the  insertion  of  a  flat-plate  w’hich 
we  had  limited  knowledge  hitherto  before.  A  strain  based  Modified  Clay  model  was 
used  to  estimate  the  associated  stress  field.  The  pore  pressure  distribution  was  com¬ 
puted  based  on  considerations  of  stress  equilibrium.  Laboratory  experiments  were 
conducted  concurrent  with  the  analytical  studies.  A  calibration  chamber  system  and 
model  flat-plate  penetrometers  were  fabricated.  Model  flat-plate  penetration  tests 
were  performed  in  laboratory  prepared  clay  specimens  using  the  chamber  system. 
Triaxial  and  consolidation  tests  were  conducted  on  the  same  clay  to  provide  reference 
properties  of  the  soil.  Field  flat-plate  penetration  test  results  from  a  Clarkson  Univer¬ 
sity  research  test  site  were  evaluated.  Based  on  the  findings,  the  following  conclusions 
are  made: 

•  A  three-dimensional  flat  plate  penetration,  as  in  cone  penetration,  results  in 
strains  and  strain  rates  much  higher  than  those  encountered  in  common  labora¬ 
tory  experiments.  Strain  rates  as  high  as  20%  per  second  are  encountered  in  the 
tip  region  of  all  penetrometers  analyzed,  as  compared  to  the  0.015%  per  second 
value  which  is  typical  to  laboratory  triaxial  testing.  The  soil  in  the  vicinity  of 
the  flat-plate  penetrometer  is  subject  to  large  strains  and  strain  reversals; 
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•  For  cylindrical  penetrometers,  ^-strain  is  the  predominant  component  as  pre¬ 
vious  studies  (Baligh,  1985)  had  indicated.  This  is,  however,  not  necessarily  the 
case  for  flat-plate  penetrometers.  The  characteristics  of  the  strain  field  during 
a  flat-plate  penetration  is  influenced  by  both  the  w/t  ratio  and  the  tip  apex 
angle.  The  combination  of  a  blunt  tip  angle  (i.  e.  ,  larger  than  20  degrees) 
and  a  w/t  of  approximately  10  or  larger  could  result  in  ^-strain  as  the  pre¬ 
dominant  strain  component  during  penetration.  This  would  render  the  lateral 
compression  theory  irrelevant  to  the  problem  it  intends  to  solve; 

•  The  basic  concept  of  linearly  relating  lateral  pressure  measurements  to  the  ex¬ 
ponential  of  plate  thickness  is  highly  questionable.  The  SPM  analyses  indicated 
that  the  predominant  strains  induced  by  the  penetration  of  a  AVStepped  Blade 
resembles  that  during  a  cyclic  direct  simple  shear  test.  The  tapered  design  is  an 
improvement  as  it  eliminates  such  cycles.  However,  a  more  serious  drawback  lies 
with  the  requirement  of  varying  w/t  ratios.  As  w/t  decreases  the  strain  paths 
deviates  further  away  from  those  of  plane  strain  penetrations,  which  is  required 
for  a  linear  exponential  relationship  between  the  lateral  pressure  measurements 
and  the  plate  thickness.  For  the  current  designs,  the  w/t  ratios  are  far  from 
being  desirable.  This  implies  that  there  is  no  theoretical  basis  for  such  a  linear 
exponential  relationship  between  pressure  and  plate  thickness,  even  if  the  soil 
is  linear  elastic.  Analytical  studies  further  indicated  that  the  pore  pressure  and 
total  normal  stress  (<7„)  peak  at  the  tip  of  the  penetrometer  regardless  of  the 
penetrometer  geometry.  The  tip  apex  angle  affects  the  peak  values  but  does 
not  alter  the  trend  of  a  sharp  decrease  of  pore  pressure  and  total  stress  as  soil 
element  passes  the  tip  of  the  penetrometer.  Computations  showed  that  the 
total  stress  remains  more  or  less  constant  as  the  soil  element  passes  the  flat- 
plate  with  increasing  thickness.  Field  stepped-blade  and  tapered  blade  tests 
confirmed  these  findings; 

•  The  flat-plate  penetration  induced  excess  pore  pressure  as  predicted  by  the 
SPM/Modified  Cam  Clay  model  system  is  significantly  lower  than  those  mea¬ 
sured  in  the  calibration  chamber.  The  plane  strain  simplification  for  the  pore 
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pressure  computation  is  obviously  a  part  of  the  reason  for  this  discrepancy.  The 
lack  of  a  more  appropriate  consideration  of  the  shear  stress  related  pore  pressure 
development  under  a  mixed  mode  of  strain  reversals  is  probably  a  more  serious 
drawback  in  the  current  approach.  Knowledge  and  available  data  base  in  this 
area  are  very  limited  at  present; 

•  For  a  simple  flat-plate  penetration  (i.e.,  the  DMT)  in  normally  consolidated 
clay,  the  magnitude  of  pore  pressure  increase  is  positively  related  to  the  soil 
rigidity  index  in  the  lateral  direction.  This  may  have  an  important  implication 
in  the  interpretation  of  DMT  tests  in  clays  to  obtain  shear  modulus  and/or 
undrained  shear  strength; 

•  The  analytical  and  experimental  studies  showed  that  the  magnitude  of  flat-plate 
penetration  induced  excess  pore  pressure  is  comparable  to  that  of  cylindrical 
penetration.  In  the  case  of  chamber  piezoblade  penetration  tests,  the  excess  pore 
pressure  diminishes  much  faster  with  distance  than  expected  for  a  cylindrical 
penetrometer  with  comparable  dimensions; 

•  The  dissipation  of  pore  pressure  at  the  center  of  a  flat-plate  is  significantly 
slower  than  expected  for  a  cylindrical  penetrometer.  The  same  has  also  been 
reported  for  field  piezoblade  tests  (Kabir  and  Lutenegger,  1990).  This  is  mostly 
due  to  a  concurrent  pore  pressure  dissipation  from  the  tip  towards  the  center 
of  the  plate.  It  is  apparent  that  placing  the  piezo  unit  within  the  tip  region  of 
flat-plate  would  make  the  interpretation  for  c/,  much  simpler; 

•  Despite  the  three-dimensional  nature  of  the  piezoblade,  the  excess  pore  pressure 
dissipation  followed  an  axisymmetric  pattern  for  the  chamber  test  performed 
in  an  anisotropically  consolidated  specimen.  Interpretation  of  the  chamber 
piezoblade  pore  pressure  dissipation  data  using  an  axisymmetric  solution  re¬ 
sulted  in  a  Ch  value  comparable  with  but  lower  than  those  obtained  from  con¬ 
solidation  tests; 

•  The  P2  readings  follow  the  same  trend  as  the  excess  pore  pressure  dissipation 
after  the  penetration  is  interrupted.  It  is  therefore  possible  to  use  P2  directly 
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from  DMT  tests  for  the  interpretation  of  q,; 

•  The  strain  path  studies  confirm  what  Hvorslev  (1949)  found  out  forty  years 
ago  that  as  far  as  disturbance  is  concerned,  it  is  the  geometry  not  the  size 
that  matters.  Any  attempt  to  improve  the  design  of  penetrometers  should 
concentrate  on  the  geometries.  Based  on  the  studies  presented  herein,  a  good 
flat-plate  penetrometer  should  have  a  w/t  exceeding  20  and  a  tip  apex  angle 
less  than  20  degrees.  The  diameter  of  the  pressure  sensor  should  be  less  than 
1  /4  of  the  plate  width.  The  sensor  lift-off  pressure  can  then  be  related  to  lateral 
compression  under  a  plane  strain  condition.  Spade  cells  (see  Figure  3.1)  used  in 
the  past  have  had  w/t  ratios  varying  from  12.7  to  25  (Massarsch,  1975;  Ladd  et 
al.  ,1979;  Tedd  and  Charles,  1981;  Fukuoka  and  Imamura,  1983).  Unfortunately, 
most  designs  did  not  take  full  advantage  of  the  relatively  favorable  w/t  ratios 
as  the  pressure  sensors  extended  the  full  width  of  the  blade.  Studies  presented 
herein  showed  that  strains  around  the  blade  corners  are  far  from  plane  strain 
and  would  undoubtedly  affect  the  pressure  measurements. 

•  If  a  flat-plate  penetrometer  (e.  g.  ,  DMT)  has  any  advantage  over  a  cylindrical 
one  (e.  g.  ,  CPT),  it  has  to  be  mechanical  or  in  operation  efficiency.  A  cone 
with  a  sharper  apex  angle  could  induce  strains  which  are  axisymmetric  and  with 
levels  comparable  to  that  of  DMT. 
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